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■ ABSTRACT 

O ■ 

Q-f We use the Chandra X-ray Observatory to study the region in the Tycho 

supernova remnant between the blast wave and the shocked ejecta interface or 
contact discontinuity. This zone contains all the history of the shock-heated gas 
and cosmic-ray acceleration in the remnant. We present for the first time evidence 
. ^ ; for significant spatial variations of the X-ray synchrotron emission in the form of 

spectral steepening from a photon index of ~ 2.6 right at the blast wave to a value 
of ~ 3.0 several arcseconds behind. We interpret this result along with the profiles 
of radio and X-ray intensity using a self-similar hydrodynamical model including 
cosmic ray backreaction that accounts for the observed ratio of radii between the 
blast wave and contact discontinuity. Two different assumptions were made about 
the post-shock magnetic field evolution: one where the magnetic field (amplified 
at the shock) is simply carried by the plasma flow and remains relatively high in 
the post-shock region [synchrotron losses limited rim case], and another where 
the amplified magnetic field is rapidly damped behind the blast wave [magnetic 
damping case]. Both cases fairly well describe the X-ray data, however both 
fail to explain the observed radio profile. The projected synchrotron emission 
leaves little room for the presence of thermal emission from the shocked ambient 
medium. This can only be explained if the pre-shock ambient medium density 
in the vicinity of the Tycho supernova remnant is below 0.6 cm -3 . 
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Subject headings: acceleration of particles — cosmic rays — ISM: individual 
(Tycho, SN 1572) — supernova remnants — X-rays: ISM 



Introduction 



One of the most remarkable discoveries made by the Chandra X-ray Observatory is that 
most of the X-ray synchrotron emission from young ejecta-dominated supernova remnants 
(SNRs) is confined to bright and geometrically thin rims at the blast wave behind which 
there is little to no evidence for thermal X-ray emission from the shocked ambient gas. 
Such thin X-ray synchrotron em itting rims are now known to be commo n in young SNRs, 
having been detecte d in Cas A ( Hughes et al. 2~OOc3: iGotthelf et al. 2001), th e Kepler SNR 



(Bambaetal 



( iBamba et al 



2003 



2005) the Tycho SNR flHwang et all 1 2002 



Warren et al. 



Long et al.ll2003l ; iRothenflug et al.l 120041 ). among others. The presence 



20051 ) and SN 1006 



of X-ray synchrotron emission requires that the SNR blast wave produces extremely high 
energy electrons (> ITeV), thereby providing strong support for efficient particle acceleration 
at high Mach number shocks. The collisionless shocks in supernova remnants are believed to 
be the main sites of production and acceleration of Galactic cosmic rays (CRs), at least up 
to the "knee" (~ 3000 TeV) of the CR spectrum. However, the new Chandra findings pose 
two important questions that have not yet been fully answered: (1) what is the underlying 
physical mechanism for the featureless thin X-ray filaments, and (2) where is the thermal 
X-ray emission from the shocked ambient medium? 

Two interpretations have been proposed to explain the thin nonthermal X-ray filaments 
seen in young SNRs: synchrotron cooling or magnetic damping. In the first case, relativistic 
electrons accelerated at the shock lose energy efficiently in the amplified post-shock magnetic 
field and after advecting or diffusing a certain distance behind the shock their synchrotron 



emission fal l s out of the X-ray b and (IBalletl 120031 ; IVink fc Laming! 120031 ; IVolk et al.l 12005 
Balletl 120061 : iParizot et al.l 120061 ). The other interpretation posits that the filaments are 



regions where the magnetic field has been amplified as well but where the rim wi dths are set 
by the damping length of the magnetic field behind the shock (jPohl et al.l 120051 ) . 



Answering the first question posed above requires that we understand the history or 
temporal evolution of both the distribution of the accelerated particles and magnetic field 
at the shock and behind. In fact, because the magnetic field strength determines the inten- 
sity of the radiative losses and subsequent changes in the energy distribution of the most 
energetic electrons, it is the most critical ingredient. It is thus necessary to understand how 
the magnetic field behaves at the shock (is it amplified?) and how it evolves downstream 
(does it decrease and, if so, at what rate?). While theoretical studies suggest that the 
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magne tic field (specifically its turbulent component) may be signific antly amplified at the 
shock (jLucek fc Bell 2000 ; Bell &: Lucek 2001 ; Vladimirov et al. 2006 ). it remains uncertain 
whether it is quic kly damped on a timescale shorter than that of the energy losses of electrons 
flPohl et alJliooih or is simply carried along by the plasma flow in the downstream region. As 
we show below these two situations lead to different predictions about the spectral variations 
of the X-ray synchrotron emission and the relative radio and X-ray synchrotron morphology 
at the blast wave. Our results will be based on a picture in which the upstream magnetic 
field is a lot stronger than is generally assumed for the interstellar medium. Assuming only 
compression (by a factor of order 4), this allows us to simply model the above mentioned 
strong magnetic field amplification. 

The other important problem is the absence of thermal X-ray emission from the shocked 
ambient medium in the region between the blast wave and the contact discontinuity in young 
ejecta-dominated SNRs. Because the blast wave also heats and compresses the ambient gas 
to X-ray emitting temperatures, a thermal X-ray component is expected in addition to the 
observed synchrotron emission. The strength of the thermal component depends on both 
the density of the ambient medium and the electronic temperature. It is unclear whether the 
lack of X-ray emission from the shocked ambient medium reflects a density or a temperature 
effect. 

The absence of a complete theory on how collisionless shocks partition their energy into 
bulk motion, thermal, and relativistic particles does not allow us to calculate the post-shock 
electronic temperature directly from dynamical quantities. If particle acceleration is very 
efficient, the shock structure is modified with respect to the case with no acceleration: the 
interaction region becomes much t hinner geometrically as well as cooler, which strongly influ- 
ences the resulting X-ray emission (IDecourchelle et al.ll2000l ; lEUison et al.ll2004l ; iDecourchelle 
20051 ). Further, even in the absence of efficient shock acceleration, the electron temperature 
can vary widely depending on whether the thermal electrons and ions at the shock front 
share the same temperature value or the same velocity distribution. In the latter case the 
electronic temperature will be much less (by the mass ratio) than the ion temperature and 
the different populations will gradually evolve toward equilibrium on the timescale set by the 
particles exchanging energy through Coulomb collisions. In fact the two cases just mentioned 
represent the extreme ranges of a continuum of possible values that we parameterize in terms 
of the electron-to-proton temperature ratio at the blast wave. There is also a spatial varia- 
tion in the X-ray intensity and spectrum of the thermal emission from the shocked ambient 
medium. Right at the shock front, the shock- heated gas has a low ionization state, which 
results in little to no K-shell line emission in the X-ray spectrum. As the gas evolves behind 
the shock, it becomes increasingly more ionized and emission lines begin to appear. Setting 
a reliable constraint on the ambient density based on the lack of the thermal X-ray emission 
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requires careful attention to the physical issues outlined above (i.e., by understanding the 
history of the thermal particles in the shocked ambient medium). 

We have chosen to address the aforementioned questions using the superb Chandra X- 
ray data of the remnant of the supernova event of 1572 that was observed by the Danish 
astronomer Tycho Brahe. Now more than 430 years later, X-ray observations of the Tycho 
SNR (hereafter, Tycho) show that the matter ejected from the explosion (the ejecta) and 
heated by the reverse shock is distributed over an almost circular shell with an angular radius 
of ~ 240". The recent advent of X-ray spectro-imagery has made it possible to identify 
the major types of nucleosynthesis products and study their spatial distribution. This, 
in turn, constrains the temperature distribution in t he shocked ejecta, the level of mixing 



and s tratification between the different ejecta layers jDecourchelle et al. 2001 ; Hwang et al. 



2003 



2006) and allows 



20021 ). the mechanism and energy of the explosion (IBadenes et al. 
us to quantify the developm ent of Raleigh- Taylor instabilities (in the X-rays, IWarren et al. 
20051 ; IVelazquez et al.lll998l . in the radio). 



In fact, the recent high resolution X-ray observations not only show a thick and clumpy 
shell of shocked ejecta but reveal the presence of thin (< 5") and smooth rims preced- 



(Hwane et al. 


2002; 


Bamba et al. 


2005; 


Warren et al. 


2005) 



sion demarcating the remnant's extent are naturally interpreted as tracing the location of 
the blast wave. Similar rims are observed in other historical SNRs such as Cas A, Kepler or 
SN 1006, but in Tycho, they are much smoother and less structured which ren ders Tycho a 



perfe ct target for the study of the rims. Their apparently featureless spectra (IHwang et al. 



20021 ) and the failure of thermal models to reproduce the observed radial profile of X-ray con- 
tinuum emission suggest that most of the X-ray r adiation from the blast wave is due to syn- 



chrotron emission by ultra-relativistic electrons (jCassam-Chena'i et al.l [2004J ; IWarren et al. 
2005l : lBalletll2006h . 



This is supported by radio observations with high an gular resolution (~ 1") t hat show a 
morphology comparable to the X-ray image as not ed by Achterberg et al. ( 1994T). althoug h 



the thin radio rims are not as highly contrasted (iDickel et al.l Il99ll ; iReynoso et al.l 119971 ) . 
By comparison Tycho has a very different appearance in the optical band. The optical 
filaments at the rim show only Ha emission from nonradiative shocks (there is no radiative 
shock emis sion from either the am bient medium or ejecta) and appear mostly in the east 
and north (jGhavamian et al.l 120001 ) . The study of the most prominent optical knot on the 



easte r n side yields a shock velo city of 1900 — 2300 km/s, independent of distance (jSmith et al. 
199ll ; iGhavamian et al.l 120011 ). The electron-to -proton temperature ra tio associated with 
this knot was constrained to be less than 0.1 (jGhavamian et al.ll200ll ). The eastern and 
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northern regions correspond to places where the SNR is possibly interacting with cold an d 
dense c louds as suggested by Hi and CO observations (jRevnoso et al.lll999l ; iLee et al.ll2004j ). 



Radio (jRevnoso et al.l 119971 ) and X-ray (IHughesI 120001 ) measurements indicate an average 
expansion rate for Tycho of 0.11 — 0.12 %yr _1 , which c orresponds to a shoc k velocity of 



3300 ± 1200 km/s at an estimated, but still debated (ISchwarz et al.lll995l). distance of 



2.3±0.8 kpc flChevalier et al.lll980l : lAlbinson et al.lll986l : ISmith et al. 



1991 



2004h . 



Ruiz-Lapuente 



I n Tycho, the q 



2005 



Warren et al 



jserv ed closeness of the contact discontinuity to the blast wave (jDecourchelle 



dio (IReynolds fc Ellison 



200 5j) and the observed concavity of the synchrotron spectrum in the ra- 



19921 ) are suggestive of efficient particle acceleration. Besides, Tycho 



being the remnant of a thermonuclear explosion, its environment is expected to be uniform 
in density and magnetic field. This uniformity is supported by the quasi-circularity and 
regularity of the rims, and by the small separation between the blast wave and the contact 
discontinuity, as opposed to Cas A where the rims appear broken all around, and far from 
the c ontact discontinuity, a result likely due to its expansion in a stellar wind (jDecourchelle 
20051 ). These characteristics all together make Tycho a particularly well-suited remnant to 
investigate the problem of the lack of thermal X-ray emission from the shocked ambient 
medium in the context of efficient particle acceleration. 

The goal of this paper is to address the above two important questions on the origin of 
the filamentary nonthermal morphology of the blast wave and lack of thermal X-ray emission 
from the shocked ambient medium. For that purpose we investigate the radial variation of 
the X-ray spectrum from the blast wave to the contact discontinuity ( §3.11 and 13. 2j) and 
compare the rim morphology in Tycho in the radio and X-ray bands ( §3.3|) . Results are 
interpreted and discussed with the use of a cosmic-ray modified hydrodynamic model of 
SNR evolution (®. 



2. Data 
2.1. X-ray 

We used the Chandra data of Tycho (obs_id 3837) observed on 2003 April 29 with the 
ACIS-I imaging spectrometer in timed exposure and faint data modes. The X-ray analysis 
was done using CIAO software (version 3.3). Standard data reduction methods were applied 
for event filtering, flare rejection, gain correction. The final exposure time amounts to 145 
ks. 

For the generation of images, spectra and PSF, we refer the reader to the on-line guide 
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of the Chandra X-ray Center websitd^ that we followed for the current analysis. Spectra 
were always adaptively grouped so th at each bin has a sign al-to- noise ratio greater than 5 a 
and we used XSPEQj (version 12.2.1. iDorman et al.ll2003l ) for the spectral modelling. 



2.2. Radio 



W e used the 1.4 GHz VLA data of Tycho observed in 1994 and 1995 (jReynoso et al. 
19971 ). These are the best currently available radio data in terms of spatial resolution (~ 1" 



Since the radio observations were done between 8 and 9 years before the X-ray obser- 
vations, we had to correct for the remnant's expansion in order to compare both the radio 
and X-ray morphologies. For radial profiles, we simply shifted the radio emission by a con- 
stant value of 3". This corresponds approximately to an expansion rate of ~ 0.35" yr _1 or 
~ 0.14% yr -1 assuming a blast wave radius of 256" . This is roughly co nsistent with the 
average expa nsion rates of 0.113% yr -1 in the radio (jReynoso et al.lll997l ) and 0.124% yr -1 
in the X-ray (jHughesI 120001 ) . Note that a shift of the radio profile by 4" would result in an 
even better match of the sharp decline of the X-ray emission in various places around the 
remnant but this corresponds to a larger expansion rate of ~ 0.18% yr^ 1 . Recently a second 
epoch Chandra observation of Tycho was approved in order to conduct a definitive X-ray 
expansion study. 



Results 



3.1. Radial variation of the synchrotron spectrum 



3.1.1. Procedure 



To search for radial spectral variations of the X-ray synchrotron emission, we selected 
regions where the thin rims at the blast wave are well defined and furthest fro m the shocked 
eject a emission. Regions with large gaps were identified using the results of IWarren et al. 
( 120051 ). while examining the 4-6 keV continuum band image (see Fig. [T]) to avoid regions 
where overlapping filaments occur at the rim. The regions are sectors chosen to have the 
largest possible extent in azimuth to include as many photons as possible. Radii were chosen 



1 See http: / /cxc. harvard.edu / ciao / threads / 
2 See http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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to best follow the local curvature of the rim. Figure [T] shows the three regions that we 
selected (W, NW and NE). Depending on the flatness of the rim, our approach leads to 
vastly different sector radii of curvature. Tabled] (see label "Rim") shows for instance that 
the resulting radius of curvature is far larger in the NE than in the W or NW, because the 
rim is nearly straight in the NE. Note that these local radii of curvature bear no relation 
to the radius of Tycho's rim, which needs to be determined with respect to the remnant's 
global center. The three azimuthally selected regions were divided radially into thin sectors 
one arcsecond wide from which X-ray spectra were extracted. 

The spectra obtained through this procedure are illustrated in Figure [2j which shows 
the evolution of the X-ray spectrum for a particular region (rim W) as one moves in radially 
from the outer boundary of the remnant toward the interior. The first few X-ray spectra 
are characterized by a well defined continuum which is sometimes associated with very 
faint K-shell emission lines of silicon and sulfur. Then, as one moves further in, these 
emission lines become more and more intense, until they clearly stand out above the X-ray 
cont inuum. While the strong emission lines can be naturally attributed to the shocked ejecta 



Decourchelle et al. 2001 



Warren et al.ll2005l ). it is not clear whether the faint lines or 



(see 

other residual emission also come from the ejecta or whether they arise from shocked ambient 
medium. 

Since we are searching for spectral variations of the X-ray synchrotron emission, we 
focus our attention on the continuum-dominated regions. There are two possible origins 
for such continuum: the thermal emission from the shocked ambient medium gas or the 
nonthermal emission from the accelerated particles. Because the morphology of the X-ray 
continuum emission and the quasi-featureless nature of the overall rim spectra cannot be 
explained by a thermal m odel in which most of the X-ray emission would come from the 
shocked ambient medium (jCassam-Chena'i et al.ll2004l ; IWarren et al.ll2005l ; iBalletl l2006l ) . we 
associate the continuum with the synchrotron emission from relativistic electrons. 



3.1.2. Blast wave and contact discontinuity 

Before proceeding to the spectral fits, we first describe how we locate the fluid discon- 
tinuities. Determining the precise radial position of the blast wave is somewhat delicate. 
Assuming spherical symmetry, the blast wave radius is given by the position where the emis- 
sion from the remnant drops to zero (after background subtraction). However, in practice 
one needs to account for the instrumental point-spread-function (PSF), which causes the 
region where the emission decreases to zero to broaden. 



- 8- 



To take account of the PSF one needs to model the emission. We have fitted a projected 
shell model with either uniform or exponential em issivity profile convo lved with the Chandra 
PSF (extracted at each particular rim), as done by lWarren et al.l (120051 ). to X-ray synchrotron 
brightness profiles (precisely, those that will be obtained in §3.1.41) . This allowed us to derive 
estimates of both the blast wave positional and width of the shell in the W, NW and NE rims 
(see Table |2]). For uniform or exponential emissivity profiles, we found an upper limit on 
the shell thickness of ~ 0.6" in the W rim with even lower values < 0.4" in the NW and NE 
rims. We note that the model profile for a thin spherical shell is less peaked than the data 
in the NW and NE rims, suggesting some deviation from our assumption of pure spherical 
symmetry. This procedure resulted in an accurate determination of the blast wave location 
for the three azimuthal regions. 



We used the results of I Warren et al.l (120051 ) to determine the ratio of blast wave to 
contact discontinuity radii and their uncertainties in the three azimuthal regions. The nu- 
merical values are R s /R c ~ 1.11 in the W rim and R s /R c — 1.09 in both rims NW and 
NE with uncertainties of order of 3 — 5%. In the numerical results we used R s /R c ~ 1.113 
(= 256"/246") for the W rim. All thes e valu es are higher than the azimuthally averaged 
value of 1.075 quoted by IWarren et al.l (120051 ). This is simply because we initially chose 
portions of the X-ray rims that were farthest from the contact discontinuity. In our analysis 
below we show how our results change as a function of R s /R c . Finally for completeness we 
note that the contact discontinuity lies some 26", 20" and 18" behind the blast wave for the 
W, NW and NE rims, respectively. In all cases the set of radial spectra we show do not 
extend this far into the remnant's interior. 



3.1.3. Power-law model 

As a first approach to search for radial spectral variations of the synchrotron emission, we 
used a phenomenological power-law to model the different rim spectra ignoring line emission. 
While this is not appropriate for the inner regions with strong emission lines (but which may 
include a certain level of synchrotron emission due to projection effects), it is a very good 
model for the quasi-featureless outer regions (see Fig. [2]). 

Figure [3] shows the best-fit parameters derived from a simple power-law model (data 
points with the symbol O) in the W, NW and NE rims. We show the power- law index (top 



3 Emission from a shock precursor could in principle alter the estimate of the blast wave radius but the 
X-ray synchrotron brightness profiles fitted with the convolved projected shell model allow little room for 
any precursor X-ray emission. 
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panel), the line-of-sight hydrogen column densitja (second panel), the X-ray brightness (from 
the fitted normalization of the power-law model at 1 keV) corrected for interstellar absorption 
(third panel) and the best-fit reduced \ 2 (bottom panel) as a function of position behind 
the blast wave. Starting from the blast wave, we see that the photon index increases over a 
few arcseconds where the X-ray synchrotron brightness is high. This is a first indication for 
radial variations of the synchrotron spectrum in the rims. Then, as we move further in, the 
photon index profile does not show a clear pattern anymore. 

In fact, we note that the photon index and the absorption appear to be highly correlated 
(see in particular the NW and NE rims). Right at the edge of the remnant where the emission 
is the brightest, the absorption is highest. This is observed in three different widely-spaced 
places around the remnant. Since there is no reason for a sudden increase of the interstellar 
absorption precisely at the position of the bright rims all around Tycho, we believe this spatial 
variation to be spurious and related to an overly simplified spectral model for the interior 
emission. We therefore fixed the absorption to a local average value for each azimuthal region 
determined by averaging the best-fit absorption values between position and position —4" 
where the rims are bright and featureless. This yields Nh = 0.73 x 10 22 cm -2 in the rim W, 
0.66 x 10 22 cm -2 in the rim NW and 0.65 x 10 22 cm" 2 in the rim NE. 

Figure [3] shows the best-fit parameters derived from a simple power-law model with a 
fixed absorption (data points marked with •) in the W, NW and NE rims. We observe an 
even more remarkable increase of the photon index (top panel) over at least ~ 5" behind 
the blast wave where the regions have a quasi-featureless spectrum. In the three rims, the 
photon index starts at ~ 2.6 — 2.8 at the shock and rises up to ~ 3.0 — 3.2. The range in 
photon index variation is about 0.3 — 0.4. The photon index profiles then reach a more or 
less uniform value further in. 

In Figure HI we demonstrate that the gradient of spectral index is robust to column 
density variations. However, it is important to have a good estimate of the interstellar ab- 
sorption because it determines the absolute value of the photon index, particularly important 
for broad-band non-thermal models. The strongest constraint other than that from rim spec- 
tra comes from the shocked ejecta X-ray emission (there is too much uncertainty associated 
with converting Hi data or optical extinction values to X-ray column densities). Therefore, 
we have fitted the ejecta spectra (see §3.1.4jl using free neutral hydrogen column density. We 
found column density values roughly consistent with those quoted above for the featureless 
rims. Again, we find a slightly larger absorption in the W (Nh = 0.72 ± 0.02 x 10 22 cm" 2 ) 



4 We used the solar abundance values from 
(WABS model in XSPEC). 



Anders fc Ebihara 



(jl982l ) for the calculation of absorption 
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than in the NW (0.66 ± 0.02 x 10 22 cm" 2 ) and NE (0.61 ± 0.02 x 10 22 cm" 2 ). Considering 
these Nh values as well as the ones derived from the simple power-law fits (see above), we 
chose to fix the column density to an approximate mid-point value of Nh = 0.7 x 10 22 cm -2 
for the subsequent analyses. 

These pure power-law fits are clearly incomplete in that they do not account for obvious 
emission lines in the spectra. The presence of these lines argues for a thermal component that 
varies across the radial sequence. Further there is the possibility that the observed steepening 
in the featureless spectra may be due to this additional thermal component, which, being 
softer than the power-law component and growing in contribution moving inward, causes 
an apparent softening of the spectral index. To test this possibility, we study two specific 
situations where the thermal emission comes either from the shocked ejecta ( §3.1.41) or from 
the shocked ambient medium ( §3.1.51) . 

3.1.4- Power-law + ejecta template 

To investigate how the observed steepening may be modified by the introduction of 
a shocked ejecta thermal component, we take for each sector (i.e., W, NW, and NE) a 
template model spectrum of the shocked ejecta and see how much of it can be included 
along with a power-law in the various radial regions. Each azimuthal region has a different 
template spectrum determined by fitting the data extracted from nearby inner regions of 
the remnant (see label "Ejecta" in Table [1]) to single- component and constant-temperature 
non-equilibrium ionization (NEI) spectral models. 

Figure [3] shows the best-fit parameters derived from a power-law model to which we add 
an ejecta component (~k data points), with the absorption held fixed, in the W, NW and NE 
rims. We can see how the profile of the photon index (top panel) is modified by comparing 
to the single power-law model fits (• data points). The spatial variations of the photon 
index profile are unchanged behind the blast wave, although there is some modification in 
the inner regions dominated by strong emission lines (see Table [3] for the rim W). As shown 
by the reduced \ 2 profile (bottom panel), the introduction of the ejecta template improves 
the fit quality in the regions where emission lines are strong as we would have expected, 
but sometimes also in regions where the emission lines are faint (see rims W and NE). This 
latter point strongly suggests that small knots of shocked ejecta have nearly reached the blast 
wave where they contribute significantly to the X-ray emission. Figures [5] and [6] illustrate 
the respective contribution of the shocked ejecta and power-law components in the NW and 
NE regions. 
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We note, however, that the radial regions where the \ 2 difference between models with 
and without the ejecta template becomes important (bottom panel of Fig. [3]), found at 
about —9" in the rims W and NW and —4" in the rim NE with respect to the blast wave, 
do not correspond to the mean position of the contact discontinuity position. At these radii 
the featureless emission from the forward shock is still the dominant broadband spectral 
component. Rather, these locations likely represent the outermost "fingers" of shocked ejecta 
resulting from the Rayleigh- Taylor instability at the contact discontinuity. On average, the 
contact discontinuity lies much further in (see §3.1.2p . 



3.1.5. Power-law + N EI model 

We consider now that the line emission could be entirely due to the shocked ambient 
gas instead of considering, as we did in the previous section, that such emission is associated 
with the ejecta material. Therefore we introduce, in addition to the power-law, a thermal 
component associated with the shocked ambient medium and see how this component affects 
the steepening of the power-law model found in our previous analysis. 

The thermal model that we choose for this sp ectral analysis is a simple NEI model with 



solar abundances (here lAnders fc Grevessd Il989l ). The parameters of this model are the 
electronic temperature kT e , the ionization age r = n e (t) dt where n e is the post-shock 
electronic density and to — t s is the flow time (i.e., the time since shock- heating), and the 
emission measure Ex oc n 2 V/ D 2 where V is the emission volume and D the distance to the 
remnant. 

When fitting the data with a power-law and the above thermal model (the absorption 
being fixed), spectral fits (not shown here) were, in the radial regions where emission lines 
are still faint, as good as or sometimes even slightly better than those obtained using the 
simple power-law model with fixed absorption. However, the model did not place the lines 
at the observed positions in the spectra of the quasi-featureless regions. In addition, there 
was no consistent pattern in the radial variation of the ionization age, while the ionization 
age of the inner regions was far too low to be consistent with the density required to explain 
the observed brightness. 

Considering this last point, we introduce a new NEI model for the shocked ambient 
medium where both the ionization age and emission measure are consistent with the same 
post-shock electronic density. We refer to this as the "self-consistent plane shock NEI model" 
or SCPNEI for short. To use this model, we need to determine the flow time and the volume 
of each emitting region, which in turn gives the ionization age and the emission measure, 
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assuming a given value for the post-shock electronic density. To estimate the flow time, we 
used a semi-analytical hydrodynamical model (to be described in more detail in §3. 2. II below) 
that is capable of reproducing the observed ratio of the blast wave to contact discontinuity 
radii. We use the same flow time over the entire volume along the line of sight at a given 
projected distance from the shock (this is of course imperfect, the true flow time is shorter 
close to the shock in the front and in the back). The volumes are determined numerically. 
Each l"-wide region has then a single-component NEI model with its own ionization age 
(fixed by the post-shock electronic density) and in which the electronic temperature is the 
only fitting parameter. 

In the following, we restrict the analysis to rim W which has the highest statistics and 
where the contact discontinuity is the farthest from the blast wave. Volumes and flow times 
at each radial position behind the blast wave (i.e., 0, —1", . . . , —15") are given in Table H] 
(see the three first columns). For each position, we looked at the variations of the photon 
index, electronic temperature and \ 2 °f our power-law plus SCPNEI model as a function 
of the post-shock electronic density. While we varied the post-shock electronic density over 
more than two orders of magnitude starting from 0.1 cm -3 , we found the power-law index 
profile to be very stable in the regions of lowest \ 2 ■ The largest differences in the spectral 
index were about 0.1 which is far too small to affect the observed steepening. 

3.1.6. Power-law + NEI model + ejecta template 

Because we showed in §3.1.41 that introducing an ejecta template improves the fit very 
significantly, the previous procedure ( §3. 1.5)) was repeated by introducing the ejecta template 
along with the power-law and SCPNEI model (see two examples in Fig. [7J. For this ultimate 
case, very similar results were obtained in terms of stability of the power-law index steepening 
(top panel of Fig. [TJ. With the introduction of the SCPNEI model, the spectral fits are 
slightly improved compared to a power-law model plus ejecta template (bottom panel of 
Fig. EJ). This is however only marginally significant (at most A% 2 = Xref ~~ X 2 = ~~ 12 for 
the inclusion of two additional parameters, where Xref i s the x 2 of the power-law plus ejecta 
template model). 

This method allows us, in addition, to constrain the electronic temperature as a function 
of the post-shock electronic density (middle panel of Fig. [7]). Because we do not convincingly 
detect the shocked ambient medium, we use the power-law plus ejecta model as a reference 
and define the allowed domain as A% 2 < (i.e., we require that the additional SCPNEI 
component not degrade the fit). As expected, when the electronic density is low (below 
~ 2cm~ 3 at position —1" and ~ 1 cm -3 at position —9"), the range of allowed temperatures 
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is not constrained but as the density gets larger, a forbidden high-temperature regime appears 
(above a few keV at position —1" and above 1 keV at position —9"). 



3.2. Where is the shocked ambient medium? 



The previous analysis has presented solid evidence for a spectral index variation in the 
synchrotron emission behind the blast wave in Tycho. We have also presented evidence for 
significant X-ray thermal emission from shocked ejecta far ahead of the contact discontinuity 
in the blast wave zone. However, evidence for a shocked ambient medium thermal component 
remains weak. The data are compatible with no such thermal component and provide 
constraints on its density and temperature. We found, in particular, that a large range 
of post-shock electronic densities were allowed when fitting the spectra at each individual 
position behind the blast wave with poor constraints on the electronic temperature for low 
densities. Even rather large density values were allowed as long as the electronic temperature 
was quite low. However there are constraints on the allowed post-shock temperatures that 
come from the hydrodynamic evolution of the remnant, as we utilize here. 



3.2.1. CR-hydro NEI model 

To better constrain the pre-shock ambient density, we use self-consistent temperature 
and hydrodynamical profiles whose parameters are determined from a semi-analytical model 
of cosmic-ray modified SNR hydrodynamics which borrows a 1-D similarity solution for the 
hydrodynamic variables. Our simulations based on this self-similar hydrodynamical calcula- 
tion are coupled with a nonlinear diffusive shock acceleration model, so th at the backreaction 



of th e particles accelerated at the blast wave is taken into account (see iDecourchelle et al. 



20001 ). This is important because these simulations, unlike those that treat the accelerated 
particles as test-particles, are able to reproduce the observed blast wave to contact discon- 
tinuity radii ratio of ~ 1.11 that we found in rim W (see §3.1.21) . This constraint results in 
an overall compression ratio, r to t, close to 6. 

For a given ambient medium density, the hydrodynamic model provides radial profiles of 
the flow time, electronic density, ionization age and mean shock temperature from the blast 
wave to the contact discontinuity (Fig. [8]). Assuming a given electron-to-proton temperature 
ratio at the blast wave, the electron temperature profile resulting from Coulomb collisional 



heating can be calculated fro m the mean shocked gas temperature and density profiles (jltoh 



19771 ; ICox &: Anderson! Il982l ). The self-consistent NEI model that we introduce and will use 



-14- 



here (hereafter the "CR-hydro NEI" model) is based on these profiles, which we divided 
into several shells to match the number of zones in rim W (see Fig. [8]). Each shell can be 
characterized by a set of average spectral parameters. 

Because of projection effects, the NEI model of a given observed zone (i.e., here a l"-wide 
sector region) is not the NEI model of the corresponding simulated shell. Indeed, each zone 
seen in projection onto the sky includes the contribution from a specific number of shells. 
Our CR-hydro NEI model in one zone is then the combination of several single-component 
constant-temperature NEI model from the shells. The volume contributions of the different 
shells to a given projected zone are computed numerically. These volumes and the post- 
shock electronic density, together with the distance (corresponding to the simulated blast 
wave radius, see Table [5]) allow us to derive the emission measure of each single NEI model. 



3.2.2. Initial parameters of the CR-hydro model 

Cosmic-ray modified hydrodynamics models were run for different values of the am- 
bient medium density (see Table [5]). There are two different ways to match the observed 
radii ratio between the blast wave and contact discontinuity in the model. We could vary 
either the injection efficiency, 77^, which is the fraction of total particle s which end up with 



suprathermal energies, o r the unshocked upstream magnetic field, B Q (lBerezhko fc Ellison 



19991 ; Ellison et al.l l2000l ). An increase in the injection efficiency increases the cosmic-ray 
pressure and then the overall compression ratio, while an increase in the upstream magnetic 
field increases the heating of the gas by the Alfv en waves in the precursor region and then 



tends to reduce the overall compression ratio (see iBerezhko fc Ellison! 1 19991 ). Since the close- 
ness of the blast wave and contact discontinuity (radii ratio of ~ 1.1) is suggestive of efficient 
particle acceleration, we choose to set the injection parameter (77^) to a high value of 1CT 3 
and adjust the upstream magnetic field for each run (see Table |5]). This case leads to lower 
post-shock gas temperatures than those derived from models with lower injection efficiency 
and hence more conservative limits on the value of the ambient density. Increasing r/ in j to 
an even higher value of 10~ 2 would in fact reduce further the post-shock pressure of the 
thermal particles albeit only by something like 15%, which would not significantly increase 
our density limits. 

We fixed the ejected mass and kinetic energy of the ejecta to 1.4 M Q and E SN = 10 51 ergs, 
respectively, which are standard values for thermonuclear SNe. The age of the SNR was fixed 
to 430 years, approximately the age of Tycho. We choose a power-law index of the initial 
power-law density profile in the ejecta, n, equal to 7 so that the expansion parameter in the 
model (m = 1 — 3/n = 4/7 ~ 0.57) is consistent with the mean expansion parameter derived 
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from the X-ray observations (0.54 ± 0.05, iHughed |2000| ). A slightly higher index n would 



produce a narrower gap between the blast wave and contact discontinuity (in the test-particle 
limit, we have for instance R s /R c = 1.14 for n = 9 compared to R s /R c = 1.18 for n = 7) 
but would clearly overestimate the remnant's expansion rate (m ~ 0.67 for n — 9) compared 
to the observations. Below, we discuss the difference between the use of a power-law and 
exponential distributions for the initial density profile in the ejecta (see §4.ip . We assume 
that the SNR evolves into an interstellar medium that is uniform in density and magnetic 
field and whose pressure is 2300 K cm" 3 . 

The self-similar model is in principle not valid as soon as the reverse shock reaches the 
core (flat density profile) of the ejecta. For n = 7 and M e j = 1.4 M Q , the mass in the ramp is 
-M e j = 0.6 M Q . The mass swept-up by the blast wave, M sw , when the reverse shock reaches 
that point is twice that for n = 7 (for no CRs), so the model is good for M sw < 1.2 M Q . 
This is reached somewhere between no = 0.1 cm -3 and 0.2 cm -3 in Table [5j This means 
that the self-similar model does not apply to higher densities. In particular, to get the same 
blast wave to contact discontinuity radii ratio when the reverse shock is inside the ejecta core 
already would presumably require a larger compression ratio to begin with. Addressing that 
goes beyond the scope of this paper and would require use of a numerical hydrodynamical 
code. 



3.2.3. Constraints on the ambient medium density 

Figure M shows the results obtained from the modelling and spectral fitting of several 
radial regions in rim W. We plot the x 2 difference between a power-law + ejecta template + 
our CR-hydro NEI model and a power-law + ejecta template as a function of the ambient 
medium density. We start from the blast wave (top panel, position 0) to the inner regions 
(bottom panel, position —9"). The different curves correspond to different initial values of 
the electron-to-proton temperature ratio at the blast wave (/3 S = 1 in solid lines and /3 S = 0.01 
in dashed lines). Note that the introduction of the CR-hydro NEI model whose parameters 
are all fixed will not necessarily always improve the fit. We only show a few regions in Figure 
|9] where the A% 2 curve goes negative below a certain density. 

The points where the \ 2 difference is null correspond to strong upper limits on the 
ambient medium density n Q . For a given electron-to-proton temperature ratio (3 S at the 
shock, these upper limits are more and more refined as we probe the innermost regions. The 
most stringent limits obtained at position —9" are no < 0.2cm" 3 for j3 s = 1 and uq < 0.3cm~ 3 
for (3 S = 0.01. We note that the case j3 s = 0.01 is not very different from the case of minimum 
equilibration (3 S = /3 min (where /3 min is the minimum electron-to-proton temperature ratio at 
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the shock given by the electron-to-proton mass ratio) for low ambient medium densities 
(below ~ 1 cm -3 ) as long as we stay near the blast wave (i.e., roughly between position 
and position —15" in Fig. |Sj). 

In Figure [101 we illustrate for the best-fit model the respective contributions of the 
power-law, shocked ejecta and shocked ambient medium for two extreme values of the 
electron-to-proton temperature ratio at the shock (/3 S = 0.01 and (3 S = 1) and an ambi- 
ent density of Uq = 0.2 cm -3 favored by the X-ray data. Table H] gives the associated photon 
index values. 

Our conclusion of this spectral analysis is that the ambient medium density must be less 
than 0.3 cm -3 . This corresponds to a lower limit on the distance to the remnant of ~ 2.8 kpc 
(Table ED which is consistent with the u pper limit range derived from optical observations 
( IChevalier et al.lll980l ; ISmith et al.lll99lT) and in good agre ement with the value derived from 
analysis of the historical light curve (jRuiz-Lapuentd 12004 ) . 



3.3. Radio and X-ray radial profiles 

3.3.1. Radio to X-ray comparison 

The first observational test to determine whether the nonthermal X-ray rims are limited 
by the magnetic field or by the energy losses of the radiating electrons consisted of searching 
for radial variations of the X-ray synchrotron spectrum and particularly for radial variations 
of its slope (see §3.ip . A second observational test, which also allows us to gain some 
insight into the spatial distribution of the magnetic field and accelerated electrons, consists 
of comparing radio and X-ray maps of the non-thermal emission at the rim (see §4. 2ft . If 
the X-ray rims are limited by the synchrotron losses, the radio synchrotron emission should 
be much broader than the X-ray synchrotron emission because the radio-emitting electrons 
are not affected by radiative losses. If the X-ray rims are magnetically limited, one expects 
radio rims as we ll, but wider and with a smaller brightness contrast between the rim and far 



behind the rim (IPohl et al.ll2005l ). 



Figure [TT] shows the radio and X-ray radial profiles of rims W, NW and NE. The X-ray 
radial profiles (data points marked with •, scale on the left) were obtained from the 2003 
Chandra image in the 4-6 keV continuum energy bands which emphasize the thin rims at the 
blast wave. The radio profiles (solid lines, scale on the right) were obtained from the 1995 
VLA image and were shifted by 3" to compensate for the remnant's expansion. The same 
spatial regions were averaged in the radio image, using the X-ray-derived radii of curvature 
(Tabled]). In this comparison, we implicitly assume that the radio flux did not change over 8 
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years. We note that the radio profiles from the NW and NE rims are similar in peak intensity, 
while the W rim is about a factor of three fainter. The comparison between the radio and 
the X-ray profiles clearly shows that the X-ray rims have a radio counterpart characterized 
by a larger width and a smaller brightness contrast between the rim and the center in the 
NW and NE rims (but not in the W rim) . This does not imply however that the X-ray rims 
are magnetically limited because the most critical constraints are the absolute radio flux and 
the ratio of X-ray to radio fluxes (see §4.2.6[) . 



3.3.2. Confidence in the radio data 



There is a contradiction between the radio profile presented in t his paper (see F ig. ITTj) . 
which were based on th e data taken i n 199 4-1995 and presented by iReynoso et all (119971 ). 
and the one obtained by IDickel et all (119911 ) in the late 1980 's (see their Fig. 2). The profile 



presented by IDickel et all (jl99ll ). a slice across the remnant that goes from S-SW to N-NE, 
shows a clear highly peaked outer filament on the N-NE side with a factor of four drop from 
this peak to the n ext valley. On the oth er hand, our attempt ^) to reproduce the profile of 
Dickel et al.l (ll99ll ) with the radio data of IReynoso et al. <jl997h shows at best a drop of only 
a factor of two. This difference is larger than what we would expect to obtain by changing 
the method to reconstruct the radio image. We did verify that pr ofiles extracted from the 
radio data we have in hand match closely the profiles published by IReynoso et al.l (119971 ). 



his differen ce suggests a possible problem in either the radio data of IDickel et al. 



(119911) or those of IReynoso et all (119971 ). In this paper, we chose to use the data presented 
by Reynoso et al. (jl997 ) because they are the most recent. In addition, in the following 
analysis, because of the above uncertainties, we will try not to use local flux values extracted 
from the radio image as inputs for models. A careful study of the relative radio and X-ray 
nonthermal emissions needs better, and better characterized, radio data. 



4. Discussion 



The two astrophysical questions that we address in this paper can be simply expressed 
as follows: (1) why is the X-ray emission so bright at the blast wave and (2) why does the 
X-ray emission fall so rapidly to faint values behind the bright rim? 



5 Note that the declination value quoted in the caption to figure 2 of 
proper sexigesimal notation and therefore may be in error. 



Dickel et al 



( 19911 ) is not given in 
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The first question aims to understand the origin of the rim morphology. Is the bright 
X-ray synchrotron emission caused by a magnetic field locally very high at the blast wave 
or by the fact that the highest energy electrons cannot travel far from their acceleration site 
without losing energy so that their emission is concentrated to very thin regions just behind 
the blast wave 

The second question aims to understand the origin of the absence of thermal X-ray 
emission from the shocked ambient medium. Is this absence caused by a low ambient medium 
density so that the thermal X-ray emission is overwhelmed by the X-ray synchrotron emission 
or by a low temperature plasma presumably resulting from efficient particle acceleration so 
that the thermal emission is hidden by interstellar absorption or shifted below the X-ray 
domain (EZD? 



4.1. The dark side of the rim 

In light of the results we present in §3.21 we conclude that the lack of thermal X- 
ray emission from the shocked ambient medium between the blast wave and the contact 
discontinuity essentially reflects a low density in the ambient medium around Tycho. 

At the blast wave, we found that the most stringent upper limit on the ambient medium 
density is about 0.9 cm -3 if electrons and protons are in temperature equilibrium and lcm -3 
if there is no equilibrium ( §3.2.31 and top panel of Fig. [9]). With these values, the thermal 
contribution from the shocked ambient medium increases far too much in the inner regions 
to be hidden by the interstellar absorption. Therefore the inner zones provide even tighter 
constraints: the ambient medium density must in fact be lower than 0.3 cm -3 . 

Our results were obtained from a spectral analysis based on different components to 
model the emission from the shock-accelerated particles, the shocked ejecta and the shocked 
ambient medium. In particular the characteristics of the shocked ambient medium were 
derived from CR-hydrodynamic models able to match the observed radii ratio between the 
blast wave and contact discontinuity as well as th e observed expansion measurements. T his 



is only possible if particle acceleration is efficient (iDecourchelld 120051 ; IWarren et al.ll2005l ). 



In that case and for any ambient density, we found from our model that the range of 
electronic temperatures allowed in the radial regions free from strong line emission is between 
5 keV and 20 keV if electron and proton temperatures are equal at the shock, and between 
0.2 keV and 2 keV if not (see Fig. [8]). This eliminates the possibility that CR hydro models, 
which satisfy the observed radii ratio, could lead to a temperature sufficiently low that the 
emission of the shocked ambient medium gets shifted to the extreme UV range. This does 
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not seem to be a viable explanation for the lack of thermal X-ray emission from the shocked 
ambient medium. 

We note that our derived upper limit value of 0.3 cm -3 , which does not depend on the 
details of the shocked ejecta emission, is inconsistent (lower by a factor ~ 3—4) with the value 
found by comparing the observed global thermal X-ray propert ies of the shocked eje cta with 



predictions from various models of thermonuclear explosions ( iBadenes et al.ll2006l ). There 
are a number of differences between this published study and the one we present here (e.g., 
the density profile of the ejecta, whether or not efficient shock acceleration is included, 
whether or not similarity solutions are used for the hydrodynamics) and both are based on 
only one-dimensional hydrodynamics. 

If we have underestimated the observed size of the gap between the blast wave and the 
contact discontinuity by a factor of two, i.e., if the ratio of radii were in reality ~1.2 (which is 
roughly the value expected in the test particle case for n = 7, then a higher ambient medium 
density would be possible since this would result in a smaller compression ratio (r tot ~ 4), 
the ionization age would increase less rapidly behind the blast wave and the volume of our 
regions would slightly decrease (because the distance would be reduced). But this would 
be no more than a factor of 2 leading to an ambient density of no ~ 0.4 — 0.6 cm~ 3 . An 
overestimate of the gap by a factor 2 (i.e., a ratio of radii of ~ 1.06 resulting in r to t ~ 10) 
would reduce the estimate on the ambient medium density by the same factor making no 
of order 0.1 — 0.2 cm" 3 . A factor of more than two error in our estimate of the size of the 
gap in either direction is unlikely. A case of lower injection efficiency (e.g., ?7i n j = 2 x 10~ 4 ) 
would lead to a higher mean shock temperature compared to the case with 77^ = 10" 3 ; this 
can only lower the density estimate (see second panel of Fig. [7]). 

Our use of a power-law initial ejecta density profile is also a possible source of un- 
certa inty. By examining density profiles generated by thermonu clear SN explosion models 
(e.g- 



Hoeflich fc Khokhlovlll996h 



Dwarkadas Chevalier! (119981 ) conclude that an exponen- 



tial profile is a reasonably good simple representation for the initial ejecta density profile. 
When evolving these profiles to the remnant stage in a uniform density environment, they 
find that the power-law and exponential profiles produce similar density and temperature 
structures in the shocked ambient medium (while they find large differences for the profiles 
in the shock heated ejecta, which we do not study here). The la rgest difference between 



the d ensity profiles in the shocked ambient medium (see Fig. 3 in iDwarkadas fc Chevalier 



19981 ) is no more than a factor of two (with the exponential case falling below the power-law 
case). Both temperature profiles are basically flat right behind the blast wave (i.e., where we 
extracted the spectra) and increase more or less rapidly only very close to the contact discon- 
tinuity (again, a region that we do not study). Although these calculations do not include 
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the effect of efficient diffuse shock acceleration, they should be indi cating roughly the level of 
differe nce between the power-law and exponential density profiles. lEllison fc Cassam-Chenai 
( 120051 ) have generated radio surface brightness profiles for thermonuclear remnants assum- 
ing power- law and exponential profiles when shock acceleration is efficient (see their Fig. 11) 
from which it is possible to estimate how the ratio R s /R c differs. Extrapolating their n = 9 
power-law results to the n = 7 power-law ejecta case we used here suggests that there is 
little difference in the ratio of radii compared to the exponential case. In summary, based on 
admittedly limited published results, we estimate that using an exponential ejecta density 
profile could result in an increase in the inferred ambient medium density of up to a factor 
of 2. 



4.2. Origin of the rim morphology 

The objective of the present discussion is to understand why most of the X-ray syn- 
chrotron emission is confined in narrow rims at the blast wave. 



4-2.1. The two interpretations 



There are currently two alternatives to explain such morphology. The first one stipulates 
that the highest energy electrons cannot travel (by advection or diffusion) far from their 
acceleration site - presumably the blast wave - without suffering from efficient energy losses 
due to synchrotron radiation so t hat their emission is concentrated to very thin regions 
(IBalletl 120031 ; IVink fc Laming 120031 ). To produce enough radiative losses this interpretation 
requires a rather high magnetic field wit hin the rims, which presuppos e s that the field mus t 
have been amplified at the blast wave (IVolk et al.l 120051 ; IBalletl 120061 ; iParizot et al.ll2006l ). 



Such turbulent amplification of the magnetic field in collisionless s 



suggested on the b asis of theoretical investigations (IBell fc Luce 



rock w aves was already 



2001) and simulations 



(jLucek fc Belli l2000l ). The second possibility recently suggested by iPohl et al.l (120051 ) is that 
the observed X-ray rims may in fact reflect the spatial distribution of the magnetic field rather 
than the spatial distribution of the high-energy electrons. This interpretation assumes also a 
certain level of amplification of the magnetic field at the blast wave (or in the precursor) and 
imposes its decrease behind, resulting from the relaxation or damping of the turbulence, on 
a timescale shorter than the characteristic time for electrons to lose energy by synchrotron 
radiation. 



The fundamental difference between these two interpretations lies in the post-shock evo- 
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lution of the magnetic field (advected or damped) or equivalently in its ability to modify the 
energy of the X-ray-emitting electrons through synchrotron losses and thereby modify their 
spatial distribution. Hence, our initial problem on the origin of the rim morphology is noth- 
ing more than a problem related to the evolution of the magnetic field behind the blast wave. 
In the following, we will attempt to determine the magnetic field characteristics/properties 
(i.e., its intensity at the shock and behind) and the parameters of the acceleration (i.e., the 
injection efficiency, the maximum energy of the shock-accelerated electrons, and the density 
ratio between the relativistic electrons and protons), by comparing the associated modelled 
synchrotron emission properties (i.e., brightness and photon index) with those derived from 
the observations. There are a large number of observational constraints which allow us to 
strongly constrain the previous parameters: the ratio of radii between the blast wave and 
contact discontinuity, the width of the X-ray rims, the X-ray spectral variations behind the 
blast wave ( §3.11) . the radio and X-ray brightness ( §3.3p . and the upper limit on the ambient 
medium density ( §3.2)) . 



4-2.2. CR-hydro model and particle spectra 

To compute the properties of the synchrotron emission, we take advantage of the CR- 
modified hydrodynamic model of SNR evolution that we used previously to estimate the 
density in the ambient medium (see 



For a given CR injection efficiency, rj^, and ambient density, no, the CR-hydro model 
provides the CR proton spectrum at the blast wave at any time. It is a piece-wise power-law 
model with an exponential cutoff at high energies: 

f p (E) = a E- r ^ exp (-£/£ p , max ) , (1) 

where a is the normalization, T is the power-law index which depends on the energy E, 
and -E Pimax is the maximum energy reached by the protons. Typically three distinct energy 



regimes with different T values are assumed (IBerezhko fe Ellison! Il999l ). The normalization 



a, is proportional to 77^ and n . The CR electron spectrum is determined by assuming a 
certain electron-to-proton density ratio at relativistic energies, K ep , which is defined as the 
ratio between the electron and proton distributions at a regime in energy wher e the protons 



are a lready relativistic but the electrons have not yet cooled radiatively (e.g., lEllison et al. 



2000). In the appropriate energy range, the CR electron spectrum is then: 

f e (E) = a K ep E~ r W exp (-E/E e>max ) , (2) 
where -E e ,max is the maximum energy reached by the electrons. K ep is left as a free parameter 
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|e.g., lEUison et al.ll2000l ; IVolk et al.ll2002l ) and will be adjusted using the X-ray data (see 
gZLi. We will obtain typically K ep ~ 10" 3 (see gXB]). 

The maximum energies E p mauX and E e mSuX contain information on the limits of the ac- 
celeration. They are set by matching either the acceleration time to the shock age or to 
the characteristic time for synchrotron losses, or by matching the upstream diffusive length 
to some fraction £ s of the shock radius (i.e., escape limitation), whichever gives the lowest 
value. We took £ s = 0.05. When the maximum energy of the electrons is limited by radiative 
losses, we have: 

3 m 2 e c 3 / (r-l)/r 1/2 1/2 

^e.max- ^3/2 V ~ + f / rB K ^2 V s, {*) 

where r and re are the overall density and magnetic field compression ratios, B2 the immedi- 
ate post-shock value of the magnetic field, V s the shock speed, and k Q = D(E e iaax ) / D B (E e)max ) 
the ratio between t he diffusion coefficient, D, and its Bohm value, Dg, both at E emax (see 



Parizot et al.ll2006l ). With typically r ~ 6 and re — 5 (see §4.2.31) . we obtain: 

£ e ,max ^ 7.3 k^ /2 B^J 2 V s>3 TeV, (4) 

where -Bioo is B 2 in units of 100 /zG, and V S:3 is V s in units of 1000 km/s. In Eq. ([2]), -Ee,max is 
left as a free parameter (via fco) and will be adjusted using the X-ray data (see §4.2.51) . We 
will obtain typically ko ~ 10 making -E e ,max ~ 10 TeV. Finally, for simplicity, we assumed 
the Bohm value and regime for the diffusion coefficient of relativistic protons. Because the 
protons affect the modelling only via their total energy density, our results are not very 
sensitive to -E^max- 

Once the particle distributions associated with their respective fluid elements are pro- 
duced, they evolve downstream experiencing adiabatic and eventually synchrotron losses (as 
described in Appendix [B]) . Then, using the electron distributions, the synchrotron emission 
within the remnant can be calculated, provided that the magnetic field structure is known 
within the remnant. 



4-2.3. Magnetic field profiles 

We present four configurations of the downstream magnetic field as a function of radius 
(normalized to the contact discontinuity) as illustrated in Figure [T2l two where the magnetic 
field is damped (left panels), and two where the magnetic field is simply advected behind 
the shock, i.e., it is passively carried by the plasma flow (right panels). Appendix IA1 explains 
how we evolve the magnetic field associated with each fluid element and then how the 
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magnetic field profiles were obtained. Anticipating the results we obtain below, for each 
magnetic field behavior (damped or advected), we show a case with an injection efficiency, 
r^inj, of 10~ 3 and pre-shock magnetic field, B , of ~ 45 uG, and one with a lower injection 
efficiency of 1.4 x 10 -4 and lower but still high pre-shock magnetic field of ~ 25 uG. Unless 
explicitly stated, these two cases are always given for an ambient density, no, of 0.2 cm -3 
and a kinetic energy of the explosion, £7sn, of 10 51 ergs. The case with 77 in j = 10~ 3 produces 
very efficient diffusive shock acceleration, and the case with ?7 in j = 1.4 x 10~ 4 yields a less 
efficient acceleration but still with some fraction of the energy flux crossing the shock going 
into relativistic particles and where the back reaction of shock-accelerated protons on the 
hydrodynamics is still important. 

The values obtained for the unshocked magnetic field, B , are clearly several times 
higher than the typical interstellar magnetic field of a few fiG. We implicitly assume that 
Bq has been already significantly amplified by some inst abilities provided, for example, by 



the s t reaming of accelerated CRs in the precursor (see iBell fc Lucekl l200ll ; iPelletier et al. 



20061 ; iMarcowith et al.l I2006T ) . Assuming the magnetic turbulence to be isotropic ahead of 
the shock, the magnetic field downstream is then larger than upstream by a factor re = 
a/(1 + 2 rf ot )/3. Since the ratio of radii between the blast wave and contact discontinuity 
in the W rim constrains the overall compression ratio, r tot , to a value of 6 (see §3.2.1j) . we 
have re — 5. This leads to an immediate post-shock magnetic field, B2, equal to ~ 215 «G 
when ?7i n j = 10~ 3 and ~ 130 uG when 77^ = 1.4 x 10~ 4 . These two sets of (?7inj, B2) will allow 
us to describe fairly well the X-ray data, i.e., the intensity of the X-ray rims (assuming a 
reasonable K ep ratio), their width and the spatial variations of the X-ray photon index. 

Figure [12] shows that, when the magnetic field is damped behind the shock (left panels), 
the final profile obtained at an age of 430 years is roughly exponentially decreasing, producing 
a magnetic filament at the blast wave. The higher the injection efficiency and immediate 
post-shock magnetic field, the sharper the filament (see Appendix lA.ip . When the magnetic 
field is advected behind the shock (right panels), the final magnetic field profile is also 
decreasing but far less so than in the damped case. These differences in the magnetic field 
profile/evolution will lead to different characteristics of the synchrotron emission. 



4-2.4- Synchrotron spectra 

Figure [13] shows the synchrotron spectrum generated at the blast wave and those pro- 
duced by several fluid elements corresponding to different observed zones (whose flow times 
are given in Table [5]) at the remnant's age of 430 years. Note that the calculated synchrotron 
emissivity, e u , was averaged over viewing angles. Our study does not include the effect of 
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magnetic field orientation. The different panels correspond to our different assumptions 
about the magnetic field evolution (damped in the left panels, advected in the right panels) 
and energy losses (only adiabatic expansion losseS in the top panels, adiabatie expansion 
plus radiative losses in the bottom panels). In those examples, the injection efficiency and 
immediate post-shock magnetic field were fixed to 1.4 x 10~ 4 and 130 /iG, and K ep to 10~ 3 
and k to 10 (left panels) or 7 (right panels). 

We illustrate the crucial role of the magnetic field (both its strength and evolution) 
on the production of the synchrotron emission. When the magnetic field is damped (left 
panels), the spectral variations of the different photon spectra come from variations in the 
strength of the magnetic field which shift the emission in frequency /energy as illustrated 
by the case with only adiabatic losses included (top- left panel). In the case of efficient 
particle acceleration considered here (77^ = 1.4 x 10~ 4 and B 2 = 130 /xG), the magnetic field 
strength diminishes rapidly as we move in from the blast wave to the interior. However, 
because the magnetic field is very large at the blast wave, its cumulative effect over time 
leads to substantial synchrotron losses and then clear changes in the synchrotron spectrum 
slope (bottom- left panel). 

On the other hand, when the magnetic field is carried by the plasma flow (right panel of 
Fig. [T3|) . differences between the case with only adiabatic losses and the one with adiabatic 
plus synchrotron losses are even more important. Because the magnetic field strength stays 
approximately constant behind the blast wave, the different spectra are almost unshifted 
in energy and therefore very similar, but in turn this produces very strong synchrotron 
losses and therefore significant spectral variations. These spectral variations are much more 
pronounced than those obtained for a decreasing magnetic field (compare the bottom panels). 



4.2.5. Methodology 

Now that the details of how we calculate the magnetic field profiles and synchrotron 
spectra are established, here we outline our methodology for relating relevant model parame- 
ters to observational constraints. The starting point is the ratio of radii, R s /R c , between the 
contact discontinuity and the blast wave, equal to 1.113 here in the W rim (see 13.1.21 . This 
provides a relation between the injection efficiency, r] in j, the post-shock magnetic field, B 2 , 
and the ambient density, n . The kinetic energy of the explosion, E SN , should be considered 
as a free parameter that also influences R s /R c . However, we freeze -Esn at 10 51 ergs, for 



6 In that case, the slope of the synchrotron spectrum reflects directly the slope of the electron spectrum 
when it was produced at the shock (see appendix |B|) . 
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simplicity. Figure [T4l shows for instance how the ratio of radii, R s /R c , varies as a function 
of B 2 for several values of rj^ at fixed n (left panel) or how it varies as a function of no for 
several values of rj in j at fixed B 2 (right panel). 

We select a value of no = 0.2 cm -3 which is fully consistent with the lack of shocked 
thermal ambient medium and derive the injection value consistent with the ratio of radii 
for a given value of magnetic field (left panel of Fig [T4l) . Then we calculate the projected 
profiles of the X-ray synchrotron brightness and spectral index, and verify that these model 
profiles are consistent with the observations. At this point, we are not making a comparison 
to the surface brightness data in flux units, but rather only to the normalized profile. The 
process is iterated using different values of B 2 in order to bracket the range of X-ray profile 
widths. We obtain reasonable fits for each magnetic field configuration (damped or advected) 
and this effectively results in a constraint on B 2 . At the same time, we can calculate the 
normalized radio profile as well as the ratio of X-ray to radio fluxes. The most important 
parameter which governs that ratio is the maximum energy reached by the electrons, £ eimax . 
We set .E'e.max (via ko) in order to get the right average X-ray spectral slope. This is done for 
each magnetic field configuration and, as we will see below, the modelled radio profiles are 
quite different. Finally we determine the electron-to-proton density ratio, K ep , by scaling 
the modelled X-ray profile to the peak values of the Chandra data at the rim. Reasonable 
values of order 10~ 3 for K cr) will be obtained. 



4-2.6. Radial profiles of the synchrotron emission 

In Figure [TBI we plot the expected line-of-sight projections of the synchrotron brightness^ 
(with both adiabatic and synchrotron losses included) in one radio (blue dotted lines) and 
one X-ray band (black solid lines) corresponding to the four magnetic field configurations 
shown in Figure [121 In addition, we show the X-ray profiles convolved with a gaussian that 
matches the Chandra PSF at 1 keV (red solid lines) which allows us a direct comparison with 
the X-ray data (points marked with •). The X-ray data points correspond to the power-law 
normalization obtained by fitting the Chandra spectra in the W rim with a power-law plus 
an ejecta template and a fixed absorption of 0.7 x 10 22 cm -2 (see §3.1.41) . and divided by the 
solid angle of each radial bin. The predicted radio profiles can be compared with the radio 
data points (marked with O in blue) which differs from those of Figure [TTJ (top panel) as 
they account for the beam size. For projection we assume the blast wave curvature in the 



7 That is (4 7r) 1 J e v (r) dl performed along the line-of-sight and where e u is the synchrotron emissivity 
per unit volume at the frequency v (in erg/s/cm _3 /Hz) which appears for instance in Figure [T51 
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out-of-sky-plane direction to be equal to the sky-plane curvature. 

When the magnetic field is damped behind the blast wave (left panels of Fig. [TBI) , the 
modelled radio (dotted lines) and X-ray (red solid lines) morphologies are very similar. Both 
radio and X-ray profiles are peaked at the blast wave, with the radio profile wider and with 
a smaller brightness contrast between the rim and far behind the rim. The limb-brightening 
of the modelled radio profiles comes closer to matching what we observe in the radio rims 
in terms of morphology (see the NW and NE rims in Fig. [TT]) . but clearly not in terms of 
intensity. When either 77^ = 10~ 3 and B 2 = 215 /iG (top-left panel) or ?7 in j = 1.4 x 1CT 4 
and B 2 = 130 /iG (bottom-left panel), the predicted radio brightness is far too low compared 
to the radio data. In the magnetic damping case, once -E e ,max (or equivalently ko) has been 
fixed in order to be consistent with the observed average X-ray slope, it is difficult to increase 
the radio intensity while keeping the same level of X-ray intensity because the ratio of X-ray 
to radio fluxes (independent of K ep ) does not depend strongly on the injection efficiency. 
However, a slight variation in k$ can increase the radio emission while keeping the same level 
of X-ray emission (provided that we adjust K ep ) and without changing the average X-ray 
slope too much. For instance, when 77 in j = 10 -3 and B 2 = 215 /iG (top-left panel), changing 
k to 14 (instead of 7) and adjusting K ep to 2.5 x 1CT 3 (instead of 2.0 x 10~ 3 ) increases the 
radio emission by a factor 1.25 (which corresponds to the ratio of K ep ) without modifying 
the X-ray intensity and only increases the averaged X-ray photon index by 0.1 (index of 2.9 
instead of 2.8). Note however that in the 1.4 x 10" 4 and B 2 = 130 /iG (bottom- 

left panel), it is not possible to make the absolute radio and X-ray fluxes at the rim and 
the average X-ray photon index all consistent with the observations just by varying k . It 
requires us to change the shape of the electron cutoff. We investigate this possibility below 
(see SJLED- 

When the magnetic field is advected behind the blast wave (right panels of Fig. ITS]) , 
the modelled radio (dotted lines) and X-ray (red solid lines) morphologies are very different. 
Contrary to the X-ray profile which is strongly peaked just behind the blast wave, the radio 
profile rises slowly to a maximum near the contact discontinuity (assuming no additional 
contribution from the ejecta). The case r/ in j = 1.4 x 10~ 4 and B 2 = 130 /iG with K ep = 
1.3 x 10~ 3 (bottom-right panel) predicts a lower radio intensity compared to the radio data 
while the case 7/ in j = 10~ 3 and B 2 = 215 /iG with K ep = 1.1 x 10~ 3 (top-right panel) predicts 
a larger radio intensity. This suggests, however, that there is an intermediate value for the 
injection efficiency that can account for the observed level of radio emission. We found that 
r/ in j = 3 x 10~ 4 with B 2 ~ 175 fiG and K ep = 1.2 x 10 -3 would provide such a good fit 
(not shown). Nevertheless, even with these best-fit parameters, we note that the observed 
rapid rise of the radio emission at the blast blast can not be reproduced accurately by the 
models. Finally, for completeness, we note that other studies based only on the integrated 



-27- 



synchrotron spectrum of Tycho from the radio to the X-ray suggest a K ep value of 4 x 10 3 , 
roughly consisten t with ours, although it was obtained with a different set of parameters 
J Yolk et al.ll2002h . 



4-2.7. Radial variations of the X-ray photon index 

In Figure [161 we plot the predicted X-ray photon index as a function of radius cor- 
responding to the four magnetic field configurations shown in Figure [121 To obtain these 
profiles, we had to construct the projected synchrotron spectrum at each radial position 
using a set of projected profiles of the synchrotron emission computed at several energies 
within the X-ray band (0.7 — 7 keV), with both adiabatic and radiative losses included. 
We show in addition the photon index profiles obtained from a set of projected brightness 
profiles at different energies convolved with a gaussian that matches the Chandra PSF at 
the appropriate energy (red solid lines). For comparison, we have plotted as data points 
the X-ray photon index obtained when fitting the rim W spectra with a power-law plus an 
ejecta template and a fixed absorption of 0.7 x 10 22 cm -2 (see §3.1.41) . Each curve in Figure 
[16] tells us how the projected synchrotron morphology changes with energy /frequency. 

Figure (TBI shows that projected X-ray photon index corresponding to the four combi- 
nations of magnetic field evolution can fairly well describe the average X-ray photon index 
measured in the W rim. A good match to the average X-ray photon index profile is obtained 
with ko = 10 in the magnetic damping case and ko = 7 in the synchrotron losses case, where 
k is defined in Eq. (J3]) ( §4.2.21) . We note that the four cases do predict different overall 
ranges in the photon index radial profile (red solid lines). The magnetic damping case (left 
panels) shows a range of ~ 0.15 when ?7 in j = 10~ 3 and B 2 = 215 /xG, and ~ 0.25 when 
7/i n j = 1.4 x 10~ 4 and B 2 = 130 /iG. The range is somewhat larger for the synchrotron losses 
limited case (right panels), where the magnetic field is advected behind the blast wave: ~ 0.20 
when ?7 in j = 10~ 3 and B 2 = 215 /iG, and ~ 0.35 when ?7 in j = 1.4 x 10~ 4 and B 2 = 130 /iG. In 
fact these ranges are roughly consistent with the range observed (~ 0.3 — 0.4), at least for 
the case ?7i n j = 1.4 x 10~ 4 and and B 2 = 130 /iG. We note however that the spatial variations 
of the modelled photon index occur over a very narrow region close to the shock, although 
the gradient of photon index is slightly broadened when the PSF is included (red solid lines). 

In fact, none of the curves is able to reproduce precisely the scale length of the observed 
spectral variations without introducing an inconsistency between the morphology of the 
projected synchrotron emission and the observed one. It can be demonstrated (in the thin 
spherical shell approximation) that most of the predicted spectral variation occurs over a 
distance corresponding to something a little less than the distance between the blast wave and 
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the location where the projected brightness falls to half its maximum value (see Appendix 
IU]) . This scale length is shorter than that over which the observed photon indices appear to 
increase at least for rim W. 



4-2.8. Dependencies on model assumptions 

In the two previous sections §4.2.61 and §4.2.71 we compared the modeled and observed 
synchrotron brightnesses and photon indices predicted by each of the the synchrotron losses 
or magnetic damping models. We found a good match to the X-ray data in terms of shape 
profiles as well as in absolute normalization. This allowed us to constrain the injection 
efficiency, rj^, the immediate post-shock magnetic field, B 2 , the ratio between the electron 
and proton distributions at relativistic energies, K ep and the maximum energy reached by the 
electrons, -E e ,max- On the other hand, for neither model did we find very good agreement to 
the radio data in terms of the absolute normalizations or shape of the rim profiles. However, 
there are several model assumptions on which our results and numerical values depend. Here 
we detail these dependencies, specifically how the limitations of the particle acceleration 
model (used to calculate the synchrotron emission properties) may impact our results, how 
important is the shape of the cutoff in the electron spectrum in our modeling, and how a 
better model for the magnetic damping can change the modeled profiles. 



There are uncertainties associated with the simple acceleration model (IBerezhko fc Ellison 



19991 ) we use to generate the particle distributions at the shock. In contrast to more exact 
Monte Carlo or kinetic shock models calculations, the simple model we employ here provides 
an analytical approximation to the particle spectrum that consists of broken power-laws with 
three slopes characterizing the low, intermediate and high energy regimes. Furthermore, the 
energy at which the slope changes at low energy is forced to be at E = m p c 2 ~ 1 GeV, which 
corresponds generally to the energy where the high energy electrons emit in the radio. More 
accurate models for the par ticle spectr um from diffusive shock acceleration (which produce 



smooth particle spectra, see lBlasill2002l ). tend to predict a relatively higher number density of 
radio-emitting electrons (by a factor of a few) around 1 GeV than at higher or lower energies 
compared to the approximate particle acceleration model used here. This will result in a 
higher radio synchrotron flux in both the magnetic damping and synchrotron losses model 
cases. However, it is difficult to quantify the amount of increase without employing these 
different calculations for the particle spectra, which is beyond the scope of this study. 

There is also uncertainty associated with the cutoff of the particle distribution function 
at high energy where the relativistic electrons emit X-rays. All previous results in this work 
were presented assuming a purely exponential cutoff in the accelerated particle spectrum 
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(see §4.2.2j) . However, dev iations from homogeneity can cause the cutoff to be narrowed 
or broadened (jPetruki 120061 ) . And in fact integrated-sp ectral fits to the sy nchrotron X-rays 
from SNRs often require that the cut-off be broadened (jEUison et al.ll200ll ). Thus we are led 
to investigate how dependent our results are to the shape of the electron cutoff. For that 
purpose, we consider a new CR electron spectrum: 



f e {E) = aK, 



cp 





E 


\ a 


_-Ee,max_ 



(5) 



where a is a number characterizing the shape of the cutoff (see lEllison et al.l 120001 ) . It 
appears that any variations in a will strongly impact the ratio of X-ray to radio fluxes as 
shown in Figure [T71 This is because most of the synchrotron radiation is emitted right 
at the blast wave, i.e., where losses by synchrotron cooling downstream (which erase any 
information on the shape of the cutoff) do not have time to modify the spectrum of the 
accelerated electrons. For instance, in the magnetic damping case when r/ in j = 1.4 x 10~ 4 
and B 2 = 130 /iG, values of a = 0.8, k = 20 and K cp = 5.5 x 10~ 3 provide the right radio 
and X-ray fluxes at the rim (bottom- left panel) and an averaged X-ray photon index of 2.9 
that is consistent with the observations. In the synchrotron losses case, we found that this 
is obtained when a = 1.3, k$ = 6 and K ep = 0.8 x 10~ 3 when 77^ = 10 -3 and B 2 = 215 /iG 
(top-right panel), and a = 0.8, k = 12 and K ep = 2.4 x 10~ 3 when ?7 in j = 1.4 x 10~ 4 and 
B 2 = 130 /iG (bottom-right panel). 

Of course, these values will change depending on the quality of the data. While there is 
little systematic uncertainty associated with the X-ray flux, the radio flux can be systemati- 
cally off by a factor of order 2 (see §3.3.2p . An underestimate (overestimate) of the observed 
radio flux would increase (decrease) the K ep values by the same factor in the models where 
a is a free parameter. This, in turn, increases (decreases) the modeled X-ray flux. To make 
it consistent with the observed X-ray flux requires modifying the parameters characterizing 
the cutoff in the electron spectrum (i.e., E e>max or k and a). An increase (decrease) of 
the observed radio flux by a factor of 2 for the same X-ray flux and average photon index 
requires that we roughly lower (increase) ko by ~ 50% and increase (lower) a by ~ 20% in 
the magnetic damping and synchrotron losses model cases (assuming rj^ = 1.4 x 10~ 4 and 
B 2 = 130 /iG). 

The point of the preceding exercise is to show that the X-ray/radio flux ratio is sensitive 
to the detailed cutoff of the particle energy spectrum, about which we have few independent 
observational constraints. While we attach little importance to the precise values of a 
(and other parameters) derived here it is comforting to note that they fall within generally 
accepted ranges. On the other hand this modest little study demonstrates that the X- 
ray/radio flux ratio by itself has little power to constrain the model parameters. 
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Finally, we consider the description of the magnetic damping model proposed by lPohl et al 



(120051 ). This model assumes a phenomenological exponential falloff in magnetic field strength 
with a characteristic length, l d , given by various possible damping mechanisms (see Appendix 
lA.ip . However, there is no reason why the dropoff should not be faster or slower. A slower 
(faster) magnetic field decay behind the shock would result in a projected profile of the 
synchrotron emission which is less (more) peaked in the radio. A slower decrease could be 
obtained if rather than the amplified magnetic field, it is for instance the amplified magnetic 
energy that is exponentially damped downstream of the shock on the spatial scale A 
further complication can be that the parallel and perpendicular components of the magnetic 
field are damped on different length scales. Taking into account those possibilities and the 
fact that the magnetic field orientation may not be negligible for the calculation of the syn- 
chrotron emissivity may change the predictions of the magnetic damping model. (Note that 
the effects due to magnetic field orientation may impact the prediction of the synchrotron 
losses model case as well). 



5. Conclusion 

The present paper addresses questions concerning the heating of the ambient gas and 
acceleration of relativistic particles at SNR blast waves. In young ejecta-dominated SNRs, 
the blast wave appears in the form of an outer geometrically thin rim where most of the 
synchrotron X-ray emission is confined. This provides strong evidence for the production 
and acceleration of cosmic-ray electrons to very high energies, right at the shock. Because 
in theory the blast wave compresses and heats the ambient gas to very high temperatures, 
a thermal X-ray component is also expected, but yet, there is little to no evidence for 
such component. The region between the blast wave and the shocked ejecta interface (or 
contact discontinuity) is actually X-ray dark. The physics associated with collisionless shocks 
in SNRs is not well understood. Where is the shocked ambient medium? What is the 
fundamental physical mechanism for the production of the thin X-ray synchrotron emitting 
rims? These are precisely the outstanding questions that we address. 

The best target for such study is probably the Tycho SNR as observed by the Chandra X- 
ray Observatory. The quasi-circularity and regularity of the X-ray synchrotron emitting rims 
in Tycho provide a very convenient framework for a combined observational and theoretical 
investigation. Our starting point was the X-ray analysis of the region between the blast wave 
and contact discontinuity, which is well resolved in Tycho. In several azimuthal regions, we 
extracted a set of spectra (between 0.7 and 7 keV) over several arcseconds as one moves in 
radially from the blast wave to the contact discontinuity. These spectra contain information 
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on the thermal and nonthermal populations that can be extracted provided that we can 
separate their respective contributions to the X-ray emission. The radial variations of the 
X-ray spectrum indicate a dominant contribution from and a softening of the synchrotron 
component behind the blast wave. We found in particular that the spectral index of the 
synchrotron component increases from a value of ~ 2.6 at the blast wave to ~ 3.0 behind 
the bright X-ray rim. The radial profiles of the X-ray synchrotron emission were compared 
to similar profiles in the radio band. These profiles rise at the blast wave in a very similar 
manner but while the X-ray emission drops rapidly the radio profile tends to remain (more 
or less) at its peak value. The radial variations of the X-ray spectrum indicate also another 
contribution (primarily Si and S line emission) from small knots of shocked ejecta that have 
nearly reached the blast wave. 

The lack of thermal contribution from the shocked ambient medium to the X-ray spec- 
trum implies, in the most general case, that the shocked ambient gas has either a low 
pre-shock density (~ 0.2 cm -3 ) with no constraint on its temperature or somewhat higher 
density but with a temperature below 1 keV. To go further, we built an emission model for 
the shocked ambient medium based on cosmic-ray hydrodynamic models which satisfy the 
observed ratio of radii (~ 1.1) between the blast wave and the contact discontinuity and 
the expansion measurements. Spectral analysis using this model indicates that the ambient 
medium density must be lower than 0.3 cm -3 , assuming a kinetic energy of the explosion of 
10 51 ergs. Systematic errors due to, for example, a more astrophysically appropriate initial 
ejecta density profile could push this limit to ~ 0.6 cm' 3 . Higher densities lead to X-ray 
emission that cannot be hidden by interstellar absorption. We found that even though the 
cosmic-ray hydrodynamic models predict that the shocked ambient gas is much less hot than 
in a pure gas shock, it is never sufficiently cool that its emission gets shifted to the extreme 
UV range. This does not seem to be a viable explanation for the lack of thermal X-ray 
emission from the shocked ambient medium. 

Much of our effort in this paper went toward modeling the intensity profiles in the radio 
and X-ray bands and the X-ray spectral variations of the synchrotron emission at the rim. 
Our goal was to determine whether the observed X-ray rims reflect the spatial distribution 
of the highest energy electrons or that of the magnetic field. The most critical ingredient in 
the modeling is the magnetic field and we have considered two scenarios for its post-shock 
evolution, assuming that it has been already amplified at the blast wave: one where the 
magnetic field is simply advected downstream from the shock and remains relatively high 
in the post-shock region, and one where the magnetic field is rapidly decreasing behind the 
shock because of the damping or relaxation of the turbulence. We refer to these as the 
synchrotron losses case and magnetic damping case, respectively. In both cases, a model 
with a cosmic-ray injection of 3 x 10" 4 , a post-shock amplification of the magnetic field up 
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to 175 fiG and an electron-to proton density at relativistic energies of order 1CT 3 accurately 
reproduces the narrow gap between the blast wave and contact discontinuity, the width and 
brightness of the X-ray synchrotron rims. This model assumes an ambient density of 0.2cm~ 3 
and a kinetic energy of the explosion of 10 51 ergs. In addition, both synchrotron losses and 
magnetic damping scenarios produce radial photon index variations that can accommodate 
the range of observed variations seen in the Chandra data. This is because, even in the 
magnetic damping case, synchrotron losses play a role in shaping the X-ray morphology 
and spectral index variations at the rim. A good match to the average X-ray photon index 
is obtained when relaxing the Bohm diffusion assumption. This constrains the diffusion 
coefficient to be ~ 7—10 times the Bohm value and implies a maximum energy of the 
electrons of 10 TeV. Right at the shock, this would correspond to a cutoff energy in the 
synchrotron spectrum of 0.3 keV. 

The grossest difference between the magnetic damping and synchrotron losses model 
cases concerns the radio synchrotron emission. Given possible systematic uncertainties as- 
sociated with the absolute radio flux, the theoretical description of the particle distribution 
functions, the model for the magnetic field decay behind the shock, the possibility of a 
broadening in the high-energy cutoff of the accelerated particle spectrum, and potentially 
the projection of the synchrotron emissivity onto the line-of-sight, we cannot reject one model 
versus the other based on the ratio of absolute X-ray and radio intensities. Nevertheless, 
we can use the shapes of the projected synchrotron emission. As regards the profile of the 
radio emission, the magnetic damping case produces a sharp rise in brightness at the blast 
wave, as observed, which the synchrotron cooling profile fails to do. One possible source 
of uncertainty in this comparison comes about because our X-ray and radio observations 
were made at widely separated times and the relative positioning of the rims in the two 
wave bands is subject to error since the remnant's angular expansion rate is still only poorly 
known. On the other hand, the synchrotron cooling profile yields a gradually rising radio 
profile behind the blast wave, a feature that is not in contradiction with the observation. 
Perhaps the actual situation is a combination of these two scenarios: a geometrically thin 
region of enhanced magnetic field right at the shock, that is only partially damped to some 
intermediate field value (i.e., 50— 100 /iG rather than the 5/xG value we assumed here) in the 
post-shock zone or where additional magnetic field has been generated by turbulent motions 
caused by the outermost pieces of ejecta (which we see in the X-ray data). This type of 
magnetic field configuration might also be more consistent with the radio profiles in the NE 
and NW rim regions (see Fig. 11). Coming up with realistic evolutionary scenarios in these 
advanced cases will require some care, but they may provide interesting constraints on the 
mysterious processes by which magnetic fields are generated at collisionless shock waves. 

As we have shown in this article, a detailed study of the X-ray and radio emissions behind 
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the SNR blast wave is crucial for understanding the nature of high Mach collisionless shocks. 
And much more remains to be done. Given the importance of the relative radio and X-ray 
rim morphologies, a key issue involves having radio and X-ray observations taken at roughly 
the same time, so that uncertainties due to the remnant's expansion can be minimized. Of 
equal importance is obtaining a reliable flux-calibrated map of the radio emission so that the 
point-to-point relationship between radio and X-ray synchrotron emission can be established 
and used to discriminate between models. The X-ray rims of Tycho, as observed by Chandra, 
are still largely unresolved. Reobservation with some portion of the rim at the prime, on-axis 
pointing location of Chandra with the narrowest PSF (rather than some 4' off-axis where 
the rims are in the current data set and the PSF is some 2") would help better determine 
their structure. Other young SNRs, observed by Chandra, are without question suitable for 
similar studies along the lines of what we have done here. Finally a deeper understanding of 
the magnetic field and synchrotron emission properties of the forward shock in Tycho would 
benefit from theoretical investigations usin g more sophisticated numerical hydrodynamical 
models of cosmic-ray modified shocks (e.g.. Ellison fc Cassam-Chena'J 120051 ) . 
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A. Evolution and profile of the magnetic field 
A.l. Damped magnetic field 



To obtain a magnetic filament as suggested by iPohl et al.l (120051 ) , we use the following 
phenomenological magnetic field profile to describe the evolution of the magnetic field in a 
given fluid element: 



B(r) = £?oc + (£ 2 - 5oo) exp 



(Al) 



where r is the position of a fluid element and r s the shock radius at the same time. B2 is the 
immediate post-shock magnetic field which later decays to with a characteristic length 
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ld- We set Boo to the value of 5 fiG as done in iPohl et al.l (120051 ). The length Id is equal to 
max {7k, Ia, k} where l^, l& and If are different damping lengths given by: 
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where Uj is the downstream flow speed (in the shock frame), ca = B 2 / V4 n p 2 the Alfven 
speeds just behind the blast wave (with B 2 and p 2 the magnetic field and mass density just 
behind the shock, respectively), t>L the turbulence velocity at the injection scale (typically 
a few hundreds of km/s) and a velocity equal to the Alfven speed for the high-/? plasma 
considered here (/3 being the pressure ratio between the gas and the magnetic field), A the 
wavelength of the turbulent magnetic field of order the Larmor radius, r L = E max /(e B 2 ), 
of the maximum energy protons and L the outer turbulence scale of order the shock radius 
r s . Since A <^ L, we have 1^ < Ia and, assuming that t>L ~ vj,, we have if ~ Ia- This yields 

—3/2 

Id = Ia oc B 2 so the higher the strength of the magnetic field just behind the shock, 
the smaller the damping length. Note that the above d amping leng t hs ha ve been slightly 
modified or rearranged compared to the formula given by IPohl et al.l (120051 1 . 



A. 2. Advected magnetic field 



We assume that the magnetic field is simply carried by the flow, frozen in the plasma, so 
that the parallel and perpendicular magnetic fie ld components, separately, ey olve conserving 
flux. Evolution equations are fully described in ICassam-Chena'i et al.l (120051 ) and references 
therein. 



B. Slope of particle spectrum after energy losses 

Let us assume that the particle differential spectrum N s can be locally described with 
a power-law with an index T s around an energy E s , at a time t s : N s = A/"(-E , s ) = K s E~ Fb . 

Due to the adiabatic expansion and radiative losses, the particle spectrum will be mod- 
ified in terms of energy (E s — > E) and density (N s — > N). We assume that this modified 
particle spectrum can be still described with a power-law with an index T, at a later time t: 
N = J\f(E) = K E' v . 
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Indepe ndently of the s hape of the particle spectrum, the change in energy and density 
is given by ([Reynolds! Il998l ): 

E = aU3 db; (B1) 

where 1/a = V/V s is the relative change in volume between time t and time t s , and a 
radiative loss term that includes both synchrotron and inverse Compton on the radiation 
field: 

0= / aB 2 cH (r)a 1/3 (r)dT (B3) 

with a = 4e 4 / (9m 4 c 7 ) a constant depending on the particle mass m, and B cS = (£> 2 + I?^.) 1 ^ 2 
an effective magnetic field which includes the magnetic field inside the remnant, B, and the 
magnetic field with energy density equal to that in the radiation field, B^. If this is the 
microwave background then B c ^ r = 3.27 pG. 

The change in slope of the particle spectrum is obtained by using Eqs (IBIj) and (1B2|) : 

T ^~T^ = {Ts ~ 2) [ 1 + ^] + 2. (B4) 

From Eq. (1B4|) . we see that adiabatic expansion only (i.e., = 0) does not cause any 
change in the slope of the particle spectrum, contrary to radiative losses. 



C. Projected photon index profile 

The projection along the line-of-sight of a radial emissivity profile S u results in a bright- 
ness profile E v of the form: 

E v (r) dz with | T e Z \ _ p2 , (CI) 

where r is the distance to the center of a sphere, p the distance between the center of the 
disk (i.e., projection of the sphere onto a plane) and the line-of-sight, and i the length of the 
line-of-sight. All quantities are expressed in units of the sphere's radius R s . 

If the e missivity decr eases from its maximum £ v>iaax with a characteristic width a v (in 
units of R a ). iBalletl (120061 ) has demonstrated that the brightness profile near the edge of the 
disk has the general form: 

B v (p) ~2R S v / 2a^,m ax g(y v ), (C2) 
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where g is a functional form and y v = (1 — p)/a v . In the case of an exponentially decreasing 
emissivity profile: 



/X 

„2 



g (x) = e~ x / e u du. (C3) 

The maximum of g occurs at ~ 0.854. 

In the same limit of small a v , the brightness toward the center is: 

B u (0)~2R s a u £ u , max . (C4) 

From Eq. (IC2I) . one can derive the slope (at given frequency z/) of the projected 
spectrum built from the projected brightness profile: 

dmB u ld\na v d\n£ VV3BX dlng(y u ) 

<W) = — » = J! "Tj C5 

amz/ 2 dmv dmv dmv 

In the exponential case, 

= ^ h(v , , with = 2 x - ^. (C6) 
amz/ 2 dmv g(x) 

Then, the slope profile a u reaches its maximum at the radius p* = 1 — a y* obtained by 
solving: 

^ =0 ^^ =0 ^i W)=0 where Hx)= _ 1+ (C7) 

dp dp a v Ayfxg{x) 

We find y* ~ 4.386. For comparison, the brightness profile decreases inwards to half its 
maximum value at y\ ~ 4.685. 

The maximum value of the photon index and its value at the edge (using Eqs IC5HC6I 
and \\mh{x) = — 1) and at the center (see Eq. IC4|) are respectively: 

x^O 

r„(„*) = i + O.W = i - ( <g± - d J^, (c») 

\ 2 J dmv dmv 

r„(i) = i + a„(i) = i_ jjEfass , (eg) 
r„(o) = i + a„(o) =i-^^- rfli ;f' max , (CIO) 

dmv dmv 

with %*) ~ 1.370. 

If a v oc \j\fv (as expected if the rims are limited by synchrotron losses), we find 
T v (p*) - r„(l) ^ 0.593 and r„(p*) - I\,(0) ~ 0.093. 
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Fig. 1. — Chandra images of Tycho. Left panel: 4-6 keV continuum band. Right panel: 
0.5-7 keV band. The image of the continuum emission emphasizes the narrow rims observed 
at the remnant's outer boundary. It is natural to associate these rims with the blast wave. 
The broadband image illustrates the closeness between the rims and the clumpy emission 
from the shocked ejecta. The regions of interest are labelled as W, NW, and NE (see Table 
[1]). Both images are corrected for exposure, vignetting and local astrophysical background 
and are displayed with a square-root scaling. 
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Fig. 2. — Spectra extracted from the rim W (see Fig. [T]) fitted with a power-law model 
and a fixed absorption of 0.7 x 10 22 cm" 2 . The numbers give the position from the blast 
wave in arcseconds. The position has been determined from an analysis detailed in §3.1.21 
and its absolute location corresponds to a sector whose inner and outer radii are 255" and 
256", respectively (see Tabled]). The dotted and dashed lines correspond to the Si xill Hea 
(1.75 - 1.94 keV) and S XV Hea (2.34 - 2.55 keV) bands, respectively. 
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Fig. 3. — Photon index T p h, X-ray absorption Nh, X-ray brightness corrected for absorption 
Bx and reduced \ 2 as a function of position behind the blast wave in the rims W (left panel), 
NW (middle panel) and NE (right panel) obtained for different spectral models: power-law 
with free absorption (O), power-law with a fixed absorption (•), power-law plus a template 
for the shocked ejecta with a fixed absorption (-k) [see Table |3] for numerical values]. The X- 
ray surface brightness profile obtained with this latter spectral model was fitted until position 
—9" with a uniform-emissivity projected shell model convolved by the Chandra PSF (best-fit 
in red line). The errors are in the range A% 2 < 2.7 (90% confidence level) on one parameter 
and are given only when % 2 /dof < 2. 
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Fig. 4. — Photon index T p h obtained by fitting the spectra of the rim W (see Fig. [2]) with a 
power-law model and for different values of a fixed absorption Nh,22 = Nh/10 22 cm -2 of: 0.5 
(O), 0.6 (•), 0.7 (□) and 0.8 (~k). The errors are in the range A% 2 < 2.7 (90% confidence 
level) on one parameter and are given only when x 2 /dof < 2. Varying the absorption by 
0.1 x 10 22 cm -2 shifts the photon index by ~ 0.15 but does not change the overall profile. In 
the rim W, the best-fit is obtained for N H ,22 — 0.7. 
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Fig. 5. — Spectra extracted from the rim NW (see Fig. [I]) fitted with a power-law model 
(dotted line) plus a template for the shocked ejecta (dashed line), the absorption being held 
fixed to 0.7 x 10 22 cm~ 2 . The sum of the two models is shown in solid line. The numbers 
give the position from the blast wave in arcseconds. 
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Fig. 6. — Same as Fig. [5] but for the rim NE (see Fig. [1]). 
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Fig. 7. — Photon index r ph (top panels), electronic temperature kT e (middle panels) and 
best-fit x 2 values (bottom panels) obtained by fitting the spectra of the rim W (see Fig. [2]) 
with a power-law model plus a template for the shocked ejecta and a NEI model where both 
the ionization age and emission measure are linked (see §3. 1.6|) as a function of the post-shock 
electronic density. This is shown for two regions behind the blast wave (left: —1" and right: 
—9"). The red lines show the values obtained with only a power- law model plus an ejecta 
template (x 2 = Xref)- The error bars plotted correspond to the range A% 2 = Xkef~X 2 < 0- 
In these models, the absorption was held fixed to 0.7 x 10 22 cm -2 . The rise in \ 2 aTj high 
density in the —9" panel is due to the lack of low temperature models in XSPEC. 




Fig. 8. — Parameters obtained from models of SNR evolution that include the backreaction 
from accelerated particles. The ratio of the blast wave to the contact discontinuity radii 
matches the one observed in the rim W (i.e., R s /R c = 1.113). Radial profiles of the flow time 
At, electronic density n e , ionization age r, mean gas temperature kT m , electronic temperature 
kT e are shown for different pre-shock ambient medium densities: no = 0.05 (solid line), 0.1 
(dotted line), 0.2 (dashed line), 0.5 (dash dot line) and 1.0 (dash dot dot line). The electronic 
temperature profile was computed assuming an initial electron-to-proton temperature ratio 
P s at the blast wave of f3 s = j3 m i n (red lines) where (3 m \ n is typically the electron-to-proton 
mass ratio (i.e., zero-equilibration) and 0.01 (black lines). The case of full temperature 
equilibration (j3 s = 1) is directly given by the mean shock temperature profile. The radial 
profiles between the blast wave and contact discontinuity were split up into a number of 
shells (vertical dotted lines) equal to the one in the observation and the average values of 
n e , r and kT e within each shell will serve as input for a thermal model (see §3.2.11) . 



-48- 



Rim W 




0.05 0.07 0.10 0.20 0.30 0.50 0.70 1.00 

Ambient Medium Density n ( cm" 3 ) 



Fig. 9. — Difference between the \ 2 obtained from a power-law + ejecta template + CR- 
hydro NEI model and the one obtained from a power-law + ejecta template (labelled as 
Xref in Table |3| as a function of the pre-shock ambient medium density n , for different 
electron-to-proton shock temperature ratios: (3 S = 1 (solid line) and 0.01 (dashed line). This 
is shown for different regions behind the blast wave (0, —1", —5", —7", —9") in the rim W. 
The filled circles (•) correspond to points where A% 2 = min A% 2 + 2.7 and the open circles 
(O) to points where A% 2 = 0. The corresponding densities are labelled for this latter case. 
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Fig. 10. — Different components of a best-fit model (to spectra extracted from the W rim) 
that includes a power-law (dotted line), a shocked ejecta template (dashed line) and a self- 
consistent NEI model (thick dash dot line), the interstellar absorption being held fixed to 
0.7 x 10 22 cm -2 . The parameters of the NEI model are derived from a CR-hydro model and 
take into account projection effects. The sum is shown in solid line. We show two cases that 
correspond to the same pre-shock ambient medium density of uq = 0.2 cm -3 but to different 
electron-to-proton temperature ratios at the shock: /3 B — 1 (red) and 0.01 (blue). 



-50- 




Fig. 11. — Comparison between radio and X-ray radial profiles of the rims W (top panel), NW 
(middle panel) and NE (bottom panel). The X-ray profiles were made from count images of 
the 4-6 keV continuum emission (• data points, scale on the left). The radio profiles (dashed 
line, scale on t he right) were obtain ed by using the 1994-1995 VLA high-resolution (~ 1") 



radio data (see lReynoso et al .11199 71 ) . Errors on the amplitude of the radio profiles are about 



10%. These profiles were shifted by 3" (solid line) to be comparable with the 2003 Chandra 
X-ray data. 
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Fig. 12. — Four post-shock magnetic field configurations that we consider to interpret the 
observed spatial and spectral variations of the synchrotron emission: two resulting from the 
damping of the turbulence (left panels) and two from the advection of the shocked plasma 
(right panels). We show the profiles for two couples of injection efficiency and immediate 
post-shock magnetic field: rj m j = 1CT 3 and B 2 = 215 /xG (top panels) and r/inj = 1.4 x 1CT 4 
and B 2 = 130 jiG (bottom panels). 
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Fig. 13. — Synchrotron spectra produced at various places within the remnant corresponding 
to the flow time At = and those indicated in Table HI The four panels correspond to 
different assumptions made on the magnetic field (MF) evolution/profile and nature of the 
energy losses. The injection efficiency, r] in j, was fixed to 1.4 x 10~ 4 and the immediate post- 
shock magnetic field, B 2 , to 130 /iG. The dashed lines indicate the 0.7-7 keV energy band. 
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Fig. 14. — Left panel: predicted ratio of radii, R s /R c , between the blast wave and the contact 
discontinuity as a function of the immediate post-shock magnetic field, B 2 , for different values 
of the injection efficiency, i]^, with a pre-shock ambient density, n , fixed to 0.2 cm -3 . The 
horizontal dotted line indicates the ratio of radii measured in the W rim of Tycho. On this 
line, the overall compression ratio is equal to 6. The two • are the values that we chose to 
illustrate in Figures [T2l [151 and [161 Above B 2 — 160 fiG, the curve with ?7i n j = 10~ 2 is above 
that with 7] in j = 10 -3 . Right panel: same but as a function of Uq with B 2 fixed to 175 /zG. 
The dashed line indicates the upper limit of 0.3 cm -3 found from the lack of thermal X-ray 
emission from the shocked ambient medium. The • indicates the lower limit on the ambient 
density (~ 0.1 cm -3 ) allowed by the CR-hydro model. Below n ~ 0.4 cm -3 , the curve with 
r^inj = 10 -2 is above that with i] in ^ = 10~ 3 . 
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Fig. 15. — Line-of sight projections of the modelled radio (1.4 GHz, blue dotted lines) and 
X-ray (1 keV, black solid lines) synchrotron brightness obtained for the post-shock magnetic 
field configurations shown in Fig. [12j The red solid lines show the X-ray profiles convolved 
with a model of the Chandra PSF. The radio (marked with O, in blue) and X-ray (marked 
with •) data points of the W rim are shown for comparison. The electron-to-proton density 
ratio at relativistic energies, K ep , was adjusted so that the model (red lines) matches the 
intensity of the X-ray rim. The radio and X-ray profiles were multiplied by 10 17 and 10 22 , 
respectively, to make them comparable on the plot. 
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Fig. 16. — X-ray photon index profiles associated with the radial distributions shown in Fig. 
[T5l We show the profiles with the Chandra PSF effects taken into account (red lines), and, 
for comparison, the Chandra data points of the W rim (marked with •). The curves are 
shown with k = 10 (left panels) and k = 7 (right panels) where k appears in Eq. ([3]) (see 
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Fig. 17. — Same as Figure fT5l but with a varying shape of the cutoff in the electron spectrum 
(i.e., where a is a free parameter in Eq. [5]). 
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Table 1: Geometrical parameters of the selected sector regions (see Fig. [T]). Angles are 
counted counterclockwise and zero starts from west. The center of each region is chosen 
in order to match the local curvature of the rims. For comparison, the nominal center of 
expansion in the X-rays is a J20 oo = 00h25ml9s, <5j 20 oo = 64°08"10" (jHughes! lioool ROSAT 
data) and the center t hat minimizes the el lipticity of the blast wave is aj2ooo = 00h25ml9.4s, 
<5j2ooo = 64°08"14.0" (IWarren et al.l 120051 ) . The last column gives the reference radial bin 
position of the blast wave (BW) that we chose for each rim. 



Region 


Region Center of Curvature 


An 


gle 


Radius 


Position 






«J2000 


<5j2000 


min 


max 


min 


max 


BW 




Eject a 










220" 


230" 




W 


Rim 


00h25ml9.5s 


64°08'14.0" 


335° 


360° 


240" 


264" 


255.5 ±0.5" 




Background 










276" 


294" 






Eject a 










178" 


188" 




NW 


Rim 


00h25ml3.0s 


64°08'11.0 ,/ 


46° 


59° 


198" 


222" 


214.5 ± 0.5" 




Background 










234" 


252" 






Eject a 










562" 


572" 




NE 


Rim 


00h24m52.5s 


64°03 , 10.0" 


111° 


119° 


582" 


606" 


599.5 ± 0.5" 




Background 










618" 


636" 
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Table 2: Best-fit blast wave (BW) radii and shell widths as obtained from different emissivity 
projected shell models convolved with a Chandra PSF. The blast wave radii are consistent 
with the radial bin that we chose as the blast wave position in Table [U 



Rim 


Uniform emissivity 


Exponential emissivity 


Fit quality 




BW Radius Width 


BW Radius 


Width 




W 


255.1 ±0.2" < 0.6" 


255.0 ± 0.2" 


< 0.6" 


Good 


NW 


214.5 ± 0.2" < 0.4" 


214.4 ± 0.2" 


< 0.1" 


Bad 


NE 


599.2 ± 0.2" < 0.4" 


599.3 ± 0.2" 


< 0.2" 


Bad 
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Table 3: Best-fit photon index T p h and \ 2 f° r different spectral models as a function of 
position behind the blast wave in the rim W: 1 - power-law with a varying absorption, 2 - 
power-law with a fixed absorption, 3 - power-law plus a template for the shocked ejecta with 
a fixed absorption. The hydrogen column density is N H = N H ,22 x 10 22 cm~ 2 . The errors are 
in the range A% 2 < 2.7 (90% confidence level) on one parameter and are given only when 
X 2 /dof < 2. 



Rim W 




1 - power- 


-law 




2 - power 


-law 


3 - power-law + ejecta template 






(free N H ) 




(N H ,22 = 


0.7) 




(N H ,22 


= 0.7) 


Position 




r p h 


X 2 (dof) 




r p h 


X 2 (dof) 




r p h 


xIef (dof) 


+ 1" 


2.89 




80 (37) 


2.57 




86 (38) 


2.52 




79 (37) 





2.67 


(2.58-2.83) 


87 (94) 


2.59 


(2.66-2.53) 


90 (95) 


2.57 


(2.51-2.64) 


91 (94) 


-1" 


2.61 


(2.52-2.69) 


180 (145) 


2.57 


(2.62-2.52) 


180 (146) 


2.55 


(2.49-2.60) 


166 (145) 


-2" 


2.70 


(2.62-2.79) 


160 (143) 


2.66 


(2.71-2.62) 


161 (144) 


2.65 


(2.60-2.70) 


150 (143) 


-3" 


2.88 


(2.78-2.99) 


164 (121) 


2.78 


(2.83-2.72) 


169 (122) 


2.77 


(2.70-2.83) 


154 (121) 


-4" 


2.64 


(2.53-2.76) 


141 (109) 


2.72 


(2.79-2.66) 


143 (110) 


2.71 


(2.64-2.78) 


131 (109) 


-5" 


2.84 


(2.73-2.97) 


108 (98) 


2.83 


(2.90-2.76) 


108 (99) 


2.82 


(2.75-2.90) 


97 (98) 


-6" 


2.86 


(2.72-3.01) 


109 (83) 


2.86 


(2.94-2.78) 


109 (84) 


2.86 


(2.75-2.96) 


101 (83) 


-7" 


2.79 


(2.65-2.93) 


82 (76) 


2.89 


(2.97-2.80) 


84 (77) 


2.89 


(2.79-2.98) 


75 (76) 


-8" 


2.97 


(2.82-3.13) 


92 (71) 


2.95 


(3.04-2.87) 


92 (72) 


2.97 


(2.84-3.08) 


69 (71) 


-9" 


2.86 


(2.72-3.01) 


116 (72) 


2.94 


(3.03-2.86) 


117 (73) 


2.96 


(2.86-3.06) 


91 (72) 


-10" 


2.77 




178 (75) 


2.81 




178 (76) 


2.77 


(2.65-2.89) 


85 (75) 


-11" 


2.88 




306 (85) 


2.92 




306 (86) 


2.90 


(2.76-3.04) 


99 (85) 


-12" 


2.79 




337 (84) 


2.79 




332 (84) 


2.66 


(2.53-2.81) 


96 (84) 


-13" 


2.95 




415 (84) 


2.89 




416 (85) 


2.82 


(2.66-3.01) 


103 (84) 


-14" 


2.88 




411 (84) 


2.94 




413 (85) 


2.93 


(2.77-3.12) 


126 (84) 


-15" 


2.93 




417 (86) 


2.93 




417 (87) 


2.93 


(2.77-3.14) 


110 (86) 
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Table 4: Best-fit photon index T p h and x 2 obtained with a model combining a power-law 
and a self-consistent non-equilibrium (NEI) model built from a CR-hydro model with solar 
abundances and n = 0.2 cm -3 , the absorption being held fixed (N H = 0.7 x 10 22 cm -2 ). The 
NEI model aims to represent the shocked ambient medium from the blast wave (position 
0) to the contact discontinuity (position —25") in the rim W. At is the flow time and the 
volume V52 is given in units of 10 52 D^ pc cm 3 where -Dkpc is the distance to the SNR in kpc. 
We show two cases with different electron-to-proton temperature ratio f3 s : 1 and 0.01. 
Rim W power-law + ejecta template + CR-hydro NEI model 

V 52 At (3 S = 1 /3 S = 0.01 



Position 




(y r ) 




Tph 




X 2 (dof) 




Tph 




X 2 (dof) 





0.6 


10 


2.59 


(2.53- 


2.67) 


84 (94) 


2.58 


(2.50- 


2.64) 


84 (94) 


-1" 


1.7 


28 


2.57 


(2.52- 


2.62) 


160 (145) 


2.55 


(2.50- 


2.60) 


159 (145) 


-2" 


2.5 


46 


2.68 


(2.63- 


2.73) 


149 (143) 


2.65 


(2.60- 


2.70) 


148 (143) 


-3" 


2.9 


64 


2.81 


(2.75- 


2.87) 


159 (121) 


2.77 


(2.71- 


2.82) 


156 (121) 


-4" 


3.4 


81 


2.77 


(2.69- 


2.84) 


123 (109) 


2.70 


(2.62- 


2.77) 


125 (109) 


-5" 


3.7 


97 


2.89 


(2.80- 


2.97) 


89 (98) 


2.81 


(2.71- 


2.88) 


89 (98) 


-6" 


4.0 


113 


2.94 


(2.84- 


3.03) 


93 (83) 


2.83 


(2.72- 


2.91) 


97 (83) 


-7" 


4.4 


130 


2.99 


(2.88- 


3.10) 


69 (76) 


2.86 


(2.75- 


2.95) 


70 (76) 


-8" 


4.6 


147 


3.09 


(2.98- 


3.22) 


52 (71) 


2.93 


(2.83- 


3.04) 


56 (71) 


-9" 


4.9 


164 


3.08 


(2.96- 


3.20) 


82 (72) 


2.91 


(2.80- 


3.01) 


81 (72) 


-10" 


5.1 


180 


2.87 


(2.76- 


3.01) 


83 (75) 


2.71 


(2.61- 


2.82) 


81 (75) 


-11" 


5.3 


196 


3.02 


(2.88- 


3.19) 


108 (85) 


2.82 


(2.70- 


2.95) 


107 (85) 


-12" 


5.5 


211 


2.75 


(2.61- 


2.92) 


97 (84) 


2.58 


(2.46- 


2.71) 


95 (84) 


-13" 


5.8 


226 


2.95 


(2.77- 


3.17) 


115 (84) 


2.71 


(2.56- 


2.86) 


110 (84) 


-14" 


6.0 


241 


3.08 


(2.90- 


3.30) 


131 (84) 


2.82 


(2.67- 


2.98) 


128 (84) 


-15" 


6.2 


256 


3.05 


(2.85- 


3.27) 


116 (86) 


2.81 


(2.67- 


2.98) 


114 (86) 
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Table 5: Immediate post-shock magnetic field B 2 , blast wave radius R s , speed V s and dis- 
tance to the remnant D obtained from a CR-hydro model for different values of the pre-shock 
ambient medium density n (assuming an injection efficiency i] in j of 1CT 3 ). We varied the 
upstream magnetic field Bq so that the blast wave to contact discontinuity radii ratio was 
consistent with the one of the observation in rim W, i.e., R s /R c = 1.113. The high value 
of the upstream magnetic field B implicitly assumes that it has been already significantly 
amplified by the CR-streaming instability. In all cases, the overall and magnetic field com- 
pression ratios are about 6 and 5, respectively. We give the distance to the SNR, D, assuming 
an angular radius of 256" for the blast wave. 



n (cm -3 ) 


B (//G) 


B 2 (//G) 


R s (pc) 


V a (km/s) 


D (kpc) 


2.00 


98 


490 


2.63 


3420 


2.12 


1.50 


89 


444 


2.74 


3570 


2.21 


1.00 


77 


384 


2.91 


3780 


2.34 


0.90 


74 


370 


2.95 


3840 


2.38 


0.80 


71 


354 


3.00 


3900 


2.42 


0.70 


68 


339 


3.06 


3980 


2.47 


0.60 


64 


320 


3.13 


4070 


2.52 


0.50 


60 


300 


3.21 


4170 


2.59 


0.40 


56 


279 


3.32 


4310 


2.67 


0.30 


50 


250 


3.45 


4490 


2.78 


0.20 


43 


215 


3.66 


4760 


2.95 


0.10 


34 


169 


4.04 


5250 


3.26 


0.09 


32 


162 


4.10 


5330 


3.31 


0.08 


31 


155 


4.17 


5420 


3.36 


0.07 


30 


149 


4.26 


5530 


3.43 


0.06 


28 


140 


4.35 


5650 


3.50 


0.05 


26 


131 


4.46 


5800 


3.60 



